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ABSTRACT 


Real-time imaging in the terahertz (THz) speetral range was achieved using an 
uncooled, 160x120 pixel infrared microbolometer camera and a milliwatt-scale quantum 
cascade laser (QCL). By replacing the camera’s original focusing optics with a 
Tsurupica-based lens and minimizing diffraction effects incurred by the QCL output 
beam, an imaging scheme was developed in which the camera’s focal plane array 
successfully detected wavelengths that are more than an order of magnitude longer than 
those for which the system is designed. Moreover, the incorporation of parabolic 
reflecting optics yielded a capability to produce high-resolution images of objects placed 
within the beam path. Despite the low laser powers employed, this scheme allows high- 
contrast imaging of various objects concealed by a wide range of nonmetallic materials— 
confirming the suitability of this technology for homeland security screening 
applications. Furthermore, the identification of relatively obscure security features in 
British currency notes suggests that Terahertz imaging could serve a future role as a 
detection mechanism against assorted counterfeiting practices. An extensive comparative 
analysis of experimental data produced using two QCLs (resonating at 2.8 and 3.6 THz) 
provides additional insight into the physics underlying these results, and suggests 
methods by which this imaging technology could be further improved. 
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EXECUTIVE SUMMARY 


Real-time imaging in the terahertz (THz) speetral range was achieved using an 
uncooled, 160x120 pixel infrared microbolometer camera and a milliwatt-scale quantum 
cascade laser (QCL). By replacing the camera’s original focusing optics with a 
Tsurupica-based lens and minimizing diffraction effects incurred by the QCL output 
beam, an imaging scheme was developed in which the camera’s focal plane array 
successfully detected wavelengths that are more than an order of magnitude longer than 
those for which the system is designed. Moreover, the incorporation of parabolic 
reflecting optics yielded a capability to produce high-resolution images of objects placed 
within the beam path. Noise equivalent temperature difference of the camera in the 
1-5 THz frequency range was calculated to be at least 3 K, confirming the need for 
external terahertz illumination when imaging in this frequency regime. After evaluating 
the effects of various operating parameters on laser performance, the QCL was found to 
perform optimally at 1.9 A in pulsed mode with a 300 kHz repetition rate and 10-20% 
duty cycle. Despite the low average laser powers employed (0.4-1.4 mW), this scheme 
allows high-contrast imaging of various objects concealed by a wide range of nonmetallic 
materials—confirming the suitability of this technology for homeland security screening 
applications. Furthermore, the identification of relatively obscure security features in 
British currency notes suggests that terahertz imaging could serve a future role as a 
detection mechanism against assorted counterfeiting practices. An extensive comparative 
analysis of experimental data produced using two QCLs (resonating at 2.8 and 3.6 THz) 
provides additional insight into the physics underlying these results, and suggests 
methods by which this imaging technology could be further improved. 

Throughout the course of this research, theoretical and experimental results were 
presented at a variety of scientific forums. The author presented various portions of this 
work at the 2007 American Physical Society Meeting, the 2007 SPIE Defense & Security 
Conference, the Ninth International Conference on Intersubband Transitions in Quantum 
Wells, SPIE Photonics West 2008, and as the keynote address of the 2008 Northern 
California Junior Sciences and Humanities Symposium. Scientific papers were published 

xiii 



in the proceedings of three of these symposia (Proc. SPIE 6893 , Proc. SPIE 6549, and 
Proc. 9'^ Int. Conf. on Intersubband Transitions in Quantum Wells), as well as in the 
March 1, 2008 issue of Optics Letters. This research was also profiled in an article 
appearing in the April 2008 issue of Laser Focus World. 
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I. 


INTRODUCTION 


A recent increase in homeland security concerns has precipitated a simultaneous 
demand for new imaging technologies. Of particular interest is the development of 
imaging systems that can quickly and efficiently detect concealed metal objects without 
posing health risks to humans. Owing to its unique spectral characteristics, radiation in 
the 0.3-10 terahertz (THz) spectral range has drawn recent attention as a new and 
potentially powerful medium for next-generation imaging technology [1]-[41]. 

Equipped with an appropriate illuminating source and sensor, terahertz imaging 
systems are capable of stand-off imaging of concealed objects and of human body 
tissue—particularly cancerous growths, which can elude x-ray based imaging detection 
[9]-[17]. Such detection agility is due to the fact that terahertz wavelengths are short 
enough to provide sub-millimeter resolution capability, yet are also long enough to 
penetrate most non-metallic materials [18]-[19]. At the same time, terahertz radiation is 
strongly reflected by nearly all metallic substances—providing a unique opportunity for 
high-contrast imaging of metal objects concealed in common materials such as fabric, 
paper, or plastic. This capability makes terahertz imaging a useful technology for 
identification of many dangerous or prohibited objects. Furthermore, many explosive 
materials absorb strongly in the terahertz spectral band—suggesting additional 
applications for stand-off spectroscopic analysis and identification of concealed 
explosives [20]-[22]. As an added benefit, the non-ionizing nature of terahertz radiation 
allows it to be used directly with humans without incurring the health hazards associated 
with other imaging technologies. For these reasons, terahertz technology holds promise 
as a multi-fold solution to what is likely the most pernicious homeland security issue 
today: fast, safe, and effective screening of passengers and their belongings at mass- 
transit accumulation points, and stand-off identification of explosive or otherwise 
hazardous materials. 

In its present state, however, terahertz-based technology falls far short of 
supporting such grand applications. Currently, most terahertz imaging systems are based 
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on either antenna-eoupled semiconductor detectors or cryogenically-cooled bolometers 
operating in the scan mode [l]-[4]. While successful, these techniques are disadvantaged 
by a lack of portability, slow response, the need for cryogenic cooling, and/or high 
expense. As a result, they are incapable of providing a complete solution to the imaging 
problem. A relatively new method of terahertz imaging involves a different approach; 
the use of microbolometer cameras. Because this technology is based upon temperature- 
driven changes in pixel resistivity that are produced by the absorption of incident photons 
(rather than electron-hole generation/recombination, as is used in most semiconductor- 
based photodetectors), the devices are not susceptible to thermal excitation and can be 
routinely operated at room temperature [42]. Furthermore, the relatively short thermal 
time constant (10 ms) of the microbolometer focal plane array (FPA) allows real-time 
imaging at television frame rates (30 Hz). Nearly all such cameras are designed and 
engineered for use in the infrared regime (typically 8-12 pm). Nevertheless, this 
dissertation reports the successful use of an optically-modified microbolometer system 
for imaging of quantum cascade laser (QCL) beams at 2.8 and 3.6 THz. Combined with 
theoretical calculations, the experimental data produced during these imaging trials 
strongly support the conclusion that this approach holds much promise as a viable 
imaging technology for future use in homeland security applications. 
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II. BACKGROUND 


Given the extent to which scientists and engineers have successfully exploited the 
electromagnetic spectrum for technological purposes, it is perhaps surprising that a region 
as broad as the 0.3-10 THz frequency range (1.0-0.3 mm wavelength range) has 
remained largely untapped to date. Indeed, as is clear from Fig. 1, this range—^which is 
often euphemistically referred to as the “terahertz gap”—is virtually the only portion of 
the spectrum between gamma rays and radio waves which has not found a wealth of 
everyday uses in contemporary society. 
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Fig. 1. Terahertz radiation (“T-rays”) with respeet to rest of the eleetromagnetic speetrum 
(from [43]). 


This void is not due to an absence of available applications in the terahertz 
regime. On the contrary, the relatively short wavelength and deep penetration properties 
of terahertz waves for most non-metals (as well as its strong reflective properties for 
metals) provide an outstanding opportunity for many types of high-resolution imaging. 
Some prominent potential applications of terahertz radiation are the visual screening of 
persons and packages for detection of concealed objects (firearms, knives, etc.), 
spectroscopic identification of explosive materials (such as lEDs in combat zones), and 
medical imaging of human tissue (e.g., for cancerous lesions). The additional fact that 
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terahertz radiation is entirely non-ionizing* and therefore relatively safe for human 
exposure further argues in favor of the utility of this frequeney regime [9]-[l 1]. 

Teehnologieal exploitation of a partieular frequeney band, however, requires the 
development of both sourees and deteetors that operate there. In the ease of terahertz 
radiation, there is a pronouneed dearth of both. This defieieney is owing to the relatively 
small photon energies assoeiated with terahertz waves, a eharaeteristie that tends to 
render eonventional methods of generation and deteetion infeasible in this regime. 
Ironieally, while terahertz radiation quanta are not suffieiently energetie to be produeed 
or deteeted using these so-ealled “optieal” deviees, they are, at the same time, excessively 
energetie for so-ealled “eleetronie” deviees that easily produee millimeter-seale and 
mierowave radiation [44]. 

Historieally, it has been the sensor, rather than the souree, that represents the 
greater engineering diffieulty in the terahertz regime [42]. Methods whieh are used with 
great sueeess for deteetion in the infrared regime (8-12 pm) are not neeessarily 
eompatible with the longer wavelengths assoeiated with terahertz radiation. By far, the 
most eommon prineiple used in infrared imaging is that of photon deteetion by the 
proeess of eleetron-hole generation/reeombination [42, [45]. In this method, photons 
whieh are passively emitted as graybody radiation from thermal sourees are individually 
absorbed by bulk semieonduetor materials; they direetly produee eleetrons and/or holes 
that ean be measured eleetronieally to determine the ineident photon flux. These 
deteetion systems ean not, however, be operated satisfaetorily at room temperature. The 
energy gap that is required for deteetion in this infrared regime is approximately 120 
meV (for k^lOpm). Although the value of kT for an objeet at room temperature is only 

26 meV, there is high risk of the eonduetion band being populated by thermal exeitation. 
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(The intrinsie earrier eoneentration for semieonduetor materials is typieally about 10 


This property derives from the faet that THz photons do not eontain suffieient energy per quanta to 
remove (i.e., ionize) an eleetron from an atom or moleeule, and therefore eannot infliet signifieant damage 
to human DNA. In spite of this faet, thermal effeets ean arise in the skin from under extended THz 
exposure, and some studies have linked development of eataraets in the human eye under prolonged 
exposure to non-ionizing radiation. Nevertheless, the health hazards posed by THz radiation are 
signifleantly less than those posed by many eonventional (e. g., x-ray) imaging teehnologies. 
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cm' at 300 K, which generates prohibitively large dark current in the detector.) 
Therefore, such detectors can perform satisfactorily only when operated at or below 
liquid nitrogen temperatures (77 K). 

Furthermore, if this same teehnology is used to engineer a deteetion device with 
energy gap eustomized to the terahertz regime, thermal effeets beeome eompletely 
intraetable under all but the lowest eryogenic temperatures. Employing a photon 
deteetion approach at 1.0 THz, for instance, would require an energy gap of only 4 meV. 
Inasmueh as this energy value is nearly an order of magnitude lower than kT at room 
temperature, it is plain that operation of sueh a device in the terahertz regime can only be 
possible by eooling the deviee to cryogenic temperatures (preferably 4 K or less, for 
whieh kT = 0.35 meV). This is not to suggest that photon detection techniques are not an 
effeetive or viable terahertz imaging methodology. Indeed, sueh eryogenie operation has 
been used quite sueeessfully for terahertz imaging. However, this sueeess eomes by 
ineurring eosts in eonvenienee, eompaetness, and eeonomies. 

By introdueing an entirely different type of teehnology, an effeetive sensing 
system ean be designed for terahertz frequencies. Mierobolometer arrays—comprised of 
miero-seale pixels with thermally-sensitive resistivities—have proven to be an effective 
basis for real-time infrared imaging systems; more reeently, researehers at the Naval 
Postgraduate Sehool (NPS) and the Massaehusetts Institute of Teehnology (MIT) have 
independently shown that they are also at least modestly responsive to radiation 
oseillating at terahertz frequeneies [34]-[41]. The goal of this dissertation researeh was to 
further investigate methods of terahertz deteetion using mierobolometer teehnology with 
quantum easeade laser (QCL) sourees, to establish teehniques for eharaeterizing and 
improving QCL output beam power, and—most importantly—to develop an optieal 
scheme under whieh various objects (particularly metallic objects) can be imaged under 
THz illumination. 
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III. THE SOURCE: QUANTUM CASCADE LASER (QCL) 


A. QCL THEORY/BACKGROUND 

As mentioned above, eurrent seientifie investigations of terahertz radiation suffer 
from a laek of available sourees. However, as shown in Fig. 2, there are still roughly a 
half-dozen methodologies that are eapable of produeing either narrowband or broadband 
terahertz radiation. While these methods are typieally bulky in both size and weight and 
not partieularly amenable to the types of applieations listed above, they eurrently provide 
the best opportunity for terahertz imaging trials available. 



Frequency (THz) 

Fig. 2. Methods of producing terahertz radiation. “DFG” refers to difference-frequency 
generation, a technique which is used to generate the THz beam studied by the author during a 
late collaboration effort with Stanford University researchers, (from [46]) 


For the research reported in this dissertation, a pair of two quantum cascade lasers 
(QCLs) tuned to individually lase at 2.8 and 3.6 THz were used for nearly all imaging 
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experiments/ The lasers were fabricated by Professor Jerome Faist’s research group at 
the University of Neuchatel in Switzerland and graciously provided, gratis, to the NPS 
Physics Department by Danielle Chamberlin and Peter Robrish of Agilent Laboratories in 
Santa Clara, CA. 


The quantum cascade laser (QCL) is a special class of semiconductor laser that 
operates upon a fundamentally different design principle than the diode lasers that are 
frequently used for visible and near-infrared light generation [47]-[48]. While both types 
of lasers use electronic transitions between energy bands as the basis for photon 
generation, the manner and physical phenomenon by which this effect is produced is 
distinctly different. Diode lasers—the traditional, more common class of semiconductor 
laser—produce laser radiation via electronic transitions between conduction and valence 
bands. Owing to the principle of lowest energy, electron population of the lower-energy 
valence band is much greater than that of the valence band; lasing is achieved in such 
devices by stimulating emission of a photon by forcing the recombination of an electron 
in the conduction band with a hole in the valence band. Because the energy of the band 
gap is uniquely fixed by selection of a particular bulk semiconductor material, each 
emitted photon has an energy equal to Eg. By electrically or optically pumping the 
semiconductor device, the normally miniscule electron and hole populations can be 
dramatically increased beyond the intrinsic carrier concentrations—allowing the 
sustainment of a population inversion, amplification by stimulated emission, and the 
output of a laser beam from the diode. Output power from the diode laser is given by the 
following equation [49]; 


P = 


rj hv (I - 
2e 



( 1 ) 


where h is Plank’s constant, I and Lh are the applied and threshold current (respectively), 
e is the electrical charge of an electron, v is the lasing frequency, and p is the total 
quantum efficiency of the device. 


t A difference-frequency generated terahertz beam was also used in the very final stages of this 
research effort; however, these imaging trials involved only the beam (rather than individual objects) and 
were primarily used only to establish the detection threshold of the detector system. 
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In contrast, as first proposed by Kazarinov and Suris in 1971 [47] and 
experimentally demonstrated by Faist et al in 1994 [48], quantum easeade lasers are 
inherently unipolar deviees. In QCLs, laser action is achieved exclusively through 
eleetron transitions between subbands within the eonduetion band. This subband 
structure is generated by fabrieating the deviee with periodie, alternating layers of 
materials with different bandgaps (e.g., AlGaAs and GaAs) and varying thiekness. The 
superlattice that is formed by this arrangement produees a one-dimensional quantum well 
confinement effeet; it is this confinement that is responsible for splitting the energies 
associated with the conduction band into discrete eleetronie subbands in whieh eleetrons 
can reside [50]. Fig. 3 illustrates the energy level strueture of a typical QCL designed to 
operate in mid-infrared wavelengths. This energy level diagram is sloped due to the 
applieation of an external electrieal bias, whieh, as will be seen, drives the passage of 
eleetrons through the deviee from one aetive region to another. 


ii 







Fig. 3. Band stmcture of a conventional three-level QCL, in which depopulation from the 
discrete lower laser level is via optical-phonon interactions (from [52]). 
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The three-level structure of the conduction band is key to the laser’s operation. 
Electrons injected into the upper subband (level 3) are dropped to level 2 within the 
active region by stimulated and spontaneous emission; each emits a photon in the 
process. To maintain population inversion between levels 3 and 2, rapid depopulation of 
level 2 is necessary. In the case of most QCLs, this is achieved by designing the energy 
difference between levels 1 and 2 to be resonant with the optical phonon energy of the 
laser material (typically GaAs). In this manner, resonant optical-phonon scattering 
ensures rapid (sub-picosecond) relaxation of electrons from level 2 to level 1 and 
therefore allows sustainment of a continual population inversion between levels 2 and 3 
at lower applied currents (i.e., reduced threshold current). 

The principal advantage of QCLs over other semiconductor lasers derives from 
the periodic design of the device and the unipolar nature of its operation. Because each 
transition is achieved without the involvement of holes, electrons are not annihilated after 
depopulating the upper laser level. This important distinction allows each electron to 
emit a single photon when transitioning from level 3 to level 2, to then rapidly transition 
from level 2 to level 1 through phonon scattering, and then tunnel from one active region 
of the quantum well structure (through the so-called injector region, identified in Fig. 3) 
and into an adjacent active region. Upon tunneling, each electron can then repeat the 
process again and again in a repetitious, cascading fashion (the process from which the 
QCL derives its name). Indeed, for a QCL fabricated with 120 periods, each electron can 
produce 120 resonant photons, thus significantly amplifying the total beam power that is 
eventually emitted from the device. 

Thus it is seen that equation 1 presented above for conventional diode lasers is, 
also, quite applicable to QCLs—save for one important distinction. The two equations 
differ only by a factor of N, the number of periods in the laser, due to the multiplicative 
effect produced by the device’s cascading structure [49]; 

p ^ N77/iv(/-4) ^2) 

2e 
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This multiplicative effect is especially important when dealing with terahertz 
lasers. Beeause the wavelengths associated with terahertz radiation are roughly an order 
of magnitude longer than those of infrared radiation, eaeh terahertz photon earries only 
one-tenth the energy of one found in the infrared regime. (This differenee is even more 
pronouneed between terahertz and visible photon energies, whieh differ by a faetor of 
roughly 500.) Moreover, the quantum effieieney faetor, p, is found to take on lower 
values at longer wavelengths due to an attendant inerease in optieal losses within the 
eavity. On the other hand, beeause materials used in QCLs tend to have wider energy 
gaps than the narrow-gap bulk erystals used in other semieonductor deviees, QCLs are 
eapable of earrying mueh larger eurrent densities than eonventional diode lasers without 
ineurring damage, t As a result of both of these factors, QCL lasers ean produce roughly 
one thousand times the output power of diode lasers operating at the same wavelength 
[49]. These advantages obviate any liabilities assoeiated with the long wavelengths of 
the terahertz regime, and make QCLs a natural and effeetive ehoiee of souree for sueh 
radiation. 

The above seheme provides an elegant and relatively effieient means of produeing 
visible and infrared radiation alike. However, beeause terahertz photon energy is less 
than the optieal phonon energies assoeiated with the device’s semieonduetor material 
(-'32 meV for GaAs), it is not suitable for use in terahertz-tuned QCLs. Sueh an energy 
arrangement would require that the gap between the first and seeond energy levels be less 
than the gap between the seeond and the third levels—in violation of quantum 
eonfinement theory and the familiar “partiele in a box” problem. As a result, QCLs 
designed for this regime (sueh as the one employed in the present study) employ a so- 
ealled “bound-to-eontinuum” transition design, in whieh the eleetrons undergo transitions 
from a bound state into a eonduetion miniband, rather than to a diserete energy state as in 
the QCL deseribed above (Fig. 4). 

i This important advantage over diode lasers should not be interpreted as an assertion that QCLs can 
tolerate arbitrarily large bias currents without damage. During the course of this research, the 2.8 and 3.6 
THz QCLs were each operated with an applied bias of roughly 500 A/cm^. However, owing to the small 
dimensions of the laser’s alloyed contact pad (200 pm x 2 mm), this current density was delivered using a 
bias of only 2.0 A. It is expected that employing currents much higher than this value will result in 
destruction of the device—^principally through irreversible Joule heating effects. 
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Fig. 4. Band structure of a bound-to-continuum quantum cascade laser, characteristic of the 
device used to produce terahertz radiation in the present research (from [53]). 


An important distinguishing feature of the bound-to-continuum design involves 
the unique nature of this downward transition. Firstly, the transition is diagonal in the 
sense that the electron must tunnel across the injection region in order to drop to the 
miniband constituting the laser’s ground state [50]. Secondly, depopulation of the ground 
state is achieved via electron thermalization rather than optical-phonon interaction. Upon 
dropping to the ground miniband and emitting terahertz photons in the process, electrons 
are rapidly removed by the external electric field, injected to the next adjacent bound 
state, and allowed to repeat the process for each of the laser’s stages. In the 3.6 THz laser 
of this study [51], the basic QCL structure was produced by molecular beam epitaxy 
(MBE) as a set of alternating Alo. 15 Gao. 85 As and (GaAs) layers of width 12.1 (3.2), 11 
(2.4), 11 (1.5), 12 (1.2), 13.8 (1.0), 16 (0.9), 16.3 (0.6), and 9 (4.5) nm. These alternating 
layers serve to produce two distinct regions within each structure: an active region (in 
which lasing occurs) and an injection region (through which electrons are resonantly 
tunneled from one active region to another). 120 periods of this structure, embedded in a 
single Plasmon waveguide with a confinement factor of r=0.33 and a waveguide loss of 
aw=4.6 cm’', produce a series of active regions with three conduction subbands that 
constitute the laser device. 
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B. TERAHERTZ QCL OPERATION 


Injection of electrons into the active region is essential to the QCL’s operation. 
Fig. 5, which illustrates the overall engineering and design of the 3.6 THz QCL device, 
provides insight into how this is accomplished. The repeating AlGaAs/GaAs quantum 
well structure (“MQW”) described above and illustrated in Fig. 3 is grown upon a 
250 pm thick substrate of semi-insulating GaAs and a 300 nm thick layer of n-doped 
GaAs [51]. The active region is approximately 14 pm thick; the ridge waveguide, 
processed by wet etching and three metallization steps, has dimensions of 2 mm by 
200 pm. Current is applied to the alloyed contact layer above the MQW region. Because 
the non-radiative lifetime of electrons in level 3 of the active region is extremely short 
(roughly 1 ps), threshold current density for the device is relatively high. However, a 
consequence of the waveguide’s relatively small dimensions is that arbitrarily high 
currents cannot be used without incurring significant Joule heating effects within the 
QCL device. 



Ge/Au/Ag/A ii 
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Fig. 5. Schematic of quantum cascade laser used in imaging experiments. The active region of 
the laser consists of 120 periods of GaAs/AlGaAs quantum wells (after [53]). 


As a result, electrical pumping of the QCL must be approached with care. 
Application of a continuous current results in rapid heating of the laser and attendant 
thermal population of the upper laser level within the conduction subband. This effect, 

which is particularly deleterious in lasers operating at terahertz frequencies (due to the 
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exceptionally small energy level separations involved), fundamentally interferes with the 
process of stimulated photon emission; the end result is a significant reduction (if not 
extinction) of QCL output power. 

To minimize thermal effects, such lasers are typically cooled to cryogenic 
temperatures (50 K or better) and operated in pulsed mode [51]. For the 3.6 THz laser 
predominantly used in the present study, a 300 kHz pulse repetition frequency (PRF) was 
typically used with 15-20% duty cycle (500-667 ns pulse width).§ It is especially 
important in such cases to maximize the efficiency of electrical power transfer by 
eliminating any impedance mismatch that exists between the pulse generation apparatus 
and the QCL. Impedance of the pulse generator and associated cables are known; QCL 
impedance, however, can only be accurately obtained by recording its I-V characteristics 
during operation. 

Matching was performed by evaluating the QCL’s dynamic impedance and using 
the data to produce a step-down transformer with an appropriate number of windings for 
use in the electrical circuit [37]. Current-voltage (I-V) characteristics of the laser were 
measured by applying electrical pulses at a PRF of 300 kHz, and 5-25% duty cycle, over 
a range of applied voltages. For each increment of applied voltage, the current flowing 
through the laser was measured with a Pearson current monitor. 

Analysis of current-voltage (I-V) data indicates that dynamic impedance of the 
QCL during operation is approximately 2.7 Q.; as seen in Fig. 6, this result is largely 
independent of duty cycle. Conversely, insofar as a larger applied voltage is required to 
produce a given current through the device at higher duty cycles, it is evident that static 
impedance increases for longer pulse widths. The reasons for this phenomenon have not 
yet been fully identified; however, it is possible to eliminate the possibility that it is 
strictly a temperature-driven effect arising from Joule heating. On the contrary. Fig. 7 
(I-V data for various cryostat temperatures) clearly demonstrates that static impedance 
actually decreases with increasing cryostat temperature for the particular QCL under 

§ It should be noted that, although optimum results were obtained using the operating parameters 
mentioned, laser action was still achievable at duty cycles as low as 5% and as high as 25% (the results of 
which will be elucidated upon in the analysis section of this dissertation). 
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investigation. This is possibly a consequence of the fact that displacement current 
through the coaxial cable capacitors increases as the duty cycle becomes smaller. 



Fig. 6. I-V data for the 3.6-TFlz QCL, operated under various duty cyeles. Dynamie 
impedanee of the deviee is given by the reeiproeal of the slope of the linear region of the curves. 



Fig. 7. I-V data of the 3.6-THz QCL, operated under 10% duty cycle, at various cryostat 
temperatures. The I-V data set labeled as “10-20 K” was collected with cold head thermostat set 
to 10 K; however, closed-cycle refrigeration was not sufficiently robust to prevent an elevation of 
temperature (by Joule heating), to as high as 20 K, at higher bias currents. 
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IV. THE SENSOR: MICROBOLOMETER EPA 


A. MICROBOLOMETER THEORY/ BACKGROUND 

Because they operate upon a fundamentally different set of physical principles, 
thermal detectors such as microbolometers enjoy a significant advantage over 
conventional infrared photon detectors with respect to the detection of infrared radiation. 
Whereas photon detectors detect radiation through the creation of electron/hole pairs 
resulting from absorbed photons (and, thus, require cooling to cryogenic temperatures to 
prevent thermal excitation across the band gap), detection in bolometers is driven solely 
by temperature changes arising from absorption of photons by the pixel membrane layer 
[42], Under normal operation in the infrared regime, graybody radiation—originating 
from objects within the camera’s field of view—falls upon the FPA; the distinct thermal 
response resulting in each pixel allows for the generation of images which are essentially 
graybody intensity maps of the scene before the camera. Since the core temperature of 
humans is typically -'10 K above that of room temperature (300K) and the NETD of most 
microbolometers is less than 0.1 K in the infrared regime [54], microbolometers can 
easily produce high-resolution images that clearly distinguish persons from their 
background environs and, indeed, from their clothing as well. Furthermore, this NETD 
value is sufficiently sensitive that high contrast is achievable in distinguishing fine 
features of the human body (e.g., the face) based upon differences in blood flow—and 
attendant temperature variations—to different regions (Fig. 8). 


17 




Fig. 8. Infrared image of the author, taken under conventional operation of the microbolometer 
camera used in this study (without optical modifications). 


Because pixel response is determined by differences in temperature, AT, of the 
pixel relative to its ambient operating temperature, the system operates effectively at 
room temperature and does not require the types of cooling infrastructure incurred by 
conventional infrared imaging systems. Indeed, it is this room-temperature operation that 
allows microbolometer cameras to achieve their characteristic properties of compactness, 
low weight, and ruggedness (Fig. 9). 



Fig. 9. The imaging system used in the present research: Infrared Solutions IR-160 infrared 
microbolometer camera (from [54]). 
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The difference between conventional bolometers and microbolometers (such as 
the type in this dissertation research) is principally one of physical scale. The underlying 
thermal physics of the two system types are fundamentally identical; however, in order to 
produce high-quality image resolution in the infrared regime, pixel pitch of a bolometer 
camera must be kept small (on the order of 50 pm). This is not possible using 
conventional engineering techniques, but is certainly achievable in the 
microelectromechanical (MEMS) domain. The fine degree of detail achievable in 
microbolometers is evident in the microphotograph of a single microbolometer pixel 
shown in Fig. 10. 



Fig. 10. The imaging system used in the present research: Infrared Solutions IR-160 infrared 
microbolometer camera (from [54]). 

The architecture that supports the upper pixel layer shown in Fig. 10 is a 
relatively simple but elegant one (Fig. 11). By using MEMS process layering of silicon 
(Si), silicon nitride (Si 3 N 4 ), vanadium oxide (VOx), and sacrificial layers of silicon 
dioxide (Si 02 ), a free-standing pixel membrane comprised of Si 3 N 4 and VOx can be 
produced that is thermally isolated from both the substrate layer (by virtue of support 
arms, which are also constructed of Si 3 N 4 ) and adjacent pixels—thus eliminating the 
threat of any lateral heat flow which would severely degrade image quality. The active 
pixel element, often referred to as the “membrane layer,” is comprised of a layer of Si 3 N 4 
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and applied VOx film. The VOx film, which exhibits a high temperature coefficient of 
resistance (TCR), is essential to device operation in the sense that it is the mechanism by 
which differences in incident radiation are detected as changes in pixel resistance. Each 
pixel is constructed such that a current bias can be applied (to allow for measuring of 
changes in voltage, AV, resulting from radiation/resistivity variations); support arms are 
designed so as to allow both electrical and thermal conduction to the substrate (which 
serves as a heat sink and electrical read-out mechanism). Typically, the FPA (and each 
pixel contained therein) is enclosed in an evacuated package, equipped with an IR- 
transmissive window, to eliminate any possible thermal conduction to the ambient 
atmosphere. Ideally, the only allowed pathway for thermal conduction from the pixel 
membrane should be through the support arms to the substrate [45]. Furthermore, if this 
exclusive conduction path is (through judicious design) minimized by designing the pixel 
with low overall thermal conductance, the principal heat loss mechanism becomes 
radiative in nature—allowing the array to operate at the background limit [45]. 



Fig. 11. Schematic of a typical MEMS-fabricated microbolometer pixel. Primary architectural 
elements (support arms and upper pixel membrane) are comprised of Si 3 N 4 . VOx is layered upon 
Si 3 N 4 membrane to provide the pixel with a high TCR response to incident radiation. Si02 is 
used as a sacrificial layer to provide spatial separation of the membrane from the substrate layer, 
thus ensuring that the only thermal conduction allowed is that from the membrane through the 
support arms (from [55]). 


20 


Operation of a microbolometer detector is driven almost entirely by the physics of 
heat flow associated with each individual pixel. Specifically, the magnitude and rate of 
temperature change within the pixel determines the sensitivity and time constant of the 
camera [45]-[55]. To properly model these thermal characteristics, the standard heat flow 
equation is applied to obtain an analytic solution for pixel temperature change, AT. 

Thermal energy incident upon the pixel in the form of external radiation (either in 
the form of passive or active illumination), P, may be dissipated through two primary 
mechanisms. Assuming convection and radiant heat losses are negligible, the thermal 
energy absorbed by the pixel partially manifests as a temperature increase in the pixel’s 
sensitive region (directly proportional to the region’s heat capacity, C, measured in J/K), 
and is partially dissipated by conductance to the support structure (proportional to the 
thermal conductance, G, measured in W/K). To the extent that the support structure 
(consisting of two micromachined arms) is attached to the constant-temperature substrate 
layer, this structure acts as a heat sink to the heated pixels [45]. 

Thus, heat flow into and out of the pixel can be represented with the first-order 
differential equation: 

;7P-C^^^-G(Ar) = 0 (3) 

dt 

where p represents the fraction of incident radiation that is actually absorbed by the 
pixel’s sensitive layer. This fraction is largely dependent upon membrane emissivity, s, 
and therefore can be significantly less than unity when the device is operated at 
wavelengths outside of the 8-14 pm design regime. Furthermore, assuming that the 
incident radiation is modulated at an angular frequency co, in general the power term will 
contain an exponential factor, exp(j(t)t). Rearranging the equation in a form that allows 
solution by application of an integration factor, one obtains: 

+ —(Ar) = ^e^"" (4) 

dt C C 

The equation can be solved analytically by first using an integration factor, u(t): 
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(5) 



where the right-hand side of the equation refleets the faet that, with inereasing time, the 
exp(-G t/C) term vanishes. Eventually, a quasi-steady state eondition is aehieved in 
whieh the pixel temperature sinusoidally inereases and deereases with the same 
frequeney at whieh the input power is modulated. When this eondition is reaehed, the 
root-mean squared temperature differenee, ATrms, is given by: 



This final elosed-form equation for ATrms makes use of the faet that the ratio 
C/G, having units of time, is in faet the thermal time eonstant of eaeh pixel (and, by 
extension, of the aggregate deteetion system as well). For purposes of pixel response 
time, then, it is elear that small values of heat eapaeitanee and large values of thermal 
eonduetanee are desirable. This is a eonsequenee of the faet that sueh a eombination of 
parameters ensure quiek dissipation of thermal energy from the pixel’s sensitive region— 
ensuring that ehanges in the ineident radiation are rapidly deteeted by the pixels. Indeed, 
it is this time eonstant whieh determines the eamera frame rate; the two parameters are 
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inversely proportional to eaeti other. A eommonly-applied rule of thumb dietates that the 
pixel time eonstant be one-third the value of the reeiproeal of the desired frame rate [45]. 
Therefore, in order to support a eamera eapable of imaging at TV-frame rates (30 Hz in 
the United States), individual pixel thermal eonstant must be no slower than 10 ms. 

However, a fast time eonstant—while desirable—is not the only (nor, even, the 
most) important metrie of mierobolometer performanee. As elearly indieated by 
equation 9, high thermal eonduetanee adversely affeets the eritieal figure of merit, AT 
(pixel temperature ehange). This is beeause high thermal eonduetanee direetly implies 
low thermal isolation. To the extent that it is through temperature ehanges produeed by 
spatial and temporal differenees in the ineident radiation that the system operates, it is in 
faet generally more important that the pixels be designed sueh that thermal eonduetivity 
is as low as possible (while still allowing a relatively short time eonstant) [45]. This is 
equivalent to produeing a pixel with exeellent thermal isolation. By thermally isolating 
the pixel sueh that both eonduetive and eonveetive transfer are minimized, radiative 
transfer beeomes the prineipal meehanism by whieh thermal energy is transferred from 
the pixel. 

To aehieve radiation-limited operation, then, eonveetive and eonduetive transfer 
to the atmosphere is minimized by evaeuating the eapsule in whieh the foeal plane array 
is paekaged. As a result, eonduetive transfer to the support strueture beeomes the 
predominant proeess of thermal dissipation from the pixel’s sensitive region—and is 
minimized through eareful seleetion of materials and dimensions used for the support 
strueture. Indeed, due to the importanee of thermal isolation to pixel performanee, design 
of a mierobolometer generally begins with the design of the support strueture. Thermal 
eonduetanee is a property that is analogous to eleetrieal eonduetanee (the reeiproeal of 
eleetrieal resistanee) in that it is direetly proportional to the eross-seetional area, but 
inversely proportional to the length, of a eonduetor. It follows that low thermal 
eonduetanee ean be obtained by fabrieating support arms of a insulator-type material 
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(Si 3 N 4 is a typical choice) with long length but small thiekness and width. Only after 
seeuring a well-isolated pixel arehiteeture are pixel responsivity and thermal response 
time eonsidered. 


With temperature eharaeteristies of the pixel, AT(Po), thus defined, it is a 
straightforward matter to develop an expression for the manner in whieh sueh 
temperature ehanges produee a pixel-speeifie output voltage to the FPA. As mentioned 
above, the pixel’s sensitive region is a thin (< 1 pm) membrane jointly eomprised of both 
an absorbing layer and a high-TCR layer. The latter layer (eommonly made of vanadium 
oxide, VOx) typieally has a temperature eoeffieient of resistanee of -.023/K [55], although 
magnitudes as high as -.045/K are attainable using alternate preparation methods [56]. 

By definition, TCR (often given the variable a) is equal to the differential ehange 
in resistanee with temperature, expressed as a pereentage [55]; 


1 dR 

a =- 

R dT 


( 10 ) 


Integrating this equation provides an expression for the total ehange in pixel resistanee 
for a given rise in temperature. As with most semieonduetors, this value is negative— 
indieating, unlike metals or supereonduetors, that resistanee deereases with inereasing 
AT. By multiplying equation 10 by JT and integrating the left- and right-hand sides of 
the resulting formula with respeet to T and R, respeetively, one obtains: 


\^=ln(R' )-ln(Ro) = ln 
i R 


f 


V 


r'" 


\adT = a(T-TJ 


To 


( 11 ) 


where the naught subseript and prime superseript denote the initial and final values, 
respeetively, of either T or R. Taking the exponential of both sides and making use of the 
faet that temperature inereases within the pixel membrane layer are suffieiently small that 
the truneated expansion exp(x)~l-i-x ean be used with aeeuraey, one obtains: 

^^Qcc(r-T„) ^ l^a(T'-TJ ( 12 ) 
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yielding, finally, an expression (with AR=R'-Ro and AT=T'-To) that is valid for relatively 
small values of AT (the regime in whieh AR ehanges linearly with AT): 


AR = a AT 


(13) 


The output voltage from a single pixel is produeed by an applied bias eurrent, ib, and the 
resistanee ehange AR that is produeed by AT in the pixel membrane. Combining Ohm’s 
law with equation 9, this signal voltage takes the form: 

= = = (14) 

G^I+co^t' 

By dividing equation 14 by Pq, then, one obtains an expression for voltage responsivity, 
3R,v—a eritieal metric for any detection system: 



^aRrj 

gJi+co't^ 


(15) 


Here, again, one sees the prime importance of a low-G value to detector performance. 
Thermal conductance of the pixel figures prominently into every important figure of 
merit that is associated with microbolometer sensitivity. Low thermal conduction is 
especially important as it relates to pixel “crosstalk.” If the pixels are not sufficiently 
spaced (or worse, are contiguous) within the focal plane array, the resulting lateral heat 
flow between pixels will result in a loss of image resolution. Indeed, it is for this reason 
that most microbolometer designs begin with a selection of dimensions and material for 
the support structure that produce the lowest G possible. After settling upon a specific 
microbolometer design—and evaluating the pixel thermal conductance that results from 
it—work begins on developing a membrane heat capacity which allows detector 
operation with a particular time constant (as determined by the relationship t=C/G). This 
heat capacity value is a function of thickness, specific heat, electrical contacts, detecting 
material, and coating of the pixel; as such, many engineering options are available for the 
purpose of securing a particular frame rates (typically 30 Hz) for the camera [45]. 

As suggested mathematically in equation 9, there are two sets of limiting cases 
that arise for ATrms —corresponding to high and low modulation frequencies. For 
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situations in which the angular modulation frequency is much faster than the camera 
frame rate (ro»l/x), AT is largely driven by the modulation frequency and heat capacity: 

(®r»l) (16) 

Gar coC 

For operating conditions in which the angular modulation frequency is much slower than 
the camera frame rate (co«l/x), AT is principally limited by the thermal conductance of 
the support structure, G. 

^ (®r«1) (17) 

As will be seen in the following sections of this dissertation, the QCL employed in this 
research is typically operated at a 300-kHz repetition rate. For a typical time constant of 
10 ms, then, it is clear that the microbolometer effectively operates within the “high- 
frequency” regime described by equation (11). Thus, equation 17 serves to identify the 
order of temperature change that can be expected to occur within the pixel in response to 
a particular incident power density. 

B. NOISE ANALYSIS 

The metric most commonly used to quantify sensitivity of a thermal imaging 
system, noise equivalent temperature difference (NETD), is defined as the temperature 
change in a large blackbody, within the field of view of an imaging system, that causes 
the signal-to-noise ratio at circuit readout to change by unity [45]. NETD is sensitive to 
both the optical characteristics of the imaging system and, especially, the pixel noise 
equivalent power (NEP). It follows that a thorough examination of the system’s noise 
sources is a necessary prerequisite for properly evaluating camera sensitivity. 

The microbolometer camera used in this research (IR-160, Infrared Solutions) 
utilizes a 160x120 pixel EPA which is sensitive to radiation falling within the 8-14 pm 
wavelength range. With over 35% of total radiative power of the 300 K blackbody 
spectrum lying within this range, the device has a relatively low noise equivalent 
temperature difference (NETD) at room temperature; this sensitivity allows for 
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generation of high-quality infrared images under passive operation (i.e., no external 
illumination is necessary). When operated in the infrared regime, the camera has a 
dynamic range of 66 dB and NETD of less than 100 mK with f/0.8 optics [54], 

These radiometric conditions, which allow for outstanding performance in the 
infrared regime, act as technical obstacles to a microbolometer detector operating in the 
terahertz regime. As graphically depicted in Fig. 12, total radiative power emitted by a 
300 K blackbody in the 1-5 THz spectral range is only one-twentieth of that associated 
with the design infrared band. Due to this stark disparity in available illumination 
power—coupled with the fact that IR-tuned microbolometer pixels exhibit sharply 
diminished absorption at terahertz frequencies—sensitivity of the FPA for the 1-5 THz 
band of interest was separately evaluated to establish whether external illumination is 
required to operate the device at such frequencies. 



Fig. 12. Planck radiation curve for a 300 K (room temperature) blackbody, showing spectral 
irradiance as a function of frequency. The area of the left and right shaded regions represents 
total power density associated with the 1-5 THz and 21.4-37.5 THz (8-14 pm) spectral band, 
respectively. 
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As a function of noise equivalent power (NEP), NEED is given by [55]; 


NETD = - f - .NEP (18) 

where Eno is the f-number of the optics used for focusing the beam onto the EPA, Adet is 
the area of a single pixel, s is the emissivity of the membrane material, and T/ is the lens 
transmissivity. E, the speetral irradianee of the souree, is given by [57]; 

HT) = (19) 

n 


where 8t is the emissivity of the illumination souree and P(T) is the exitance of the souree 
in the terahertz frequeney range of interest. Elsing Planek’s radiation law and the faet 
that, at 300 K, hv< kT for the 1-5 THz range [57]; 


P{T) « I 


2nkT 


dv' 


27tkT j 
3c' 


( 20 ) 


where Vc is the high frequency cutoff of the terahertz region (5 THz). This formula yields 
a total incident power density over the 1-5 THz region of approximately 12 W/m at 
300 K—signifieantly less than that associated with the 8-14 pm wavelength range 
(-170 W/m ). Thus, the NETD of the deteetor in the terahertz regime can be suceinctly 
expressed as; 

NETD = - ^—jNEP (21) 

^de,^^T,kE,V^ 

The microbolometer eamera used in the present study employs an f/1 lens with 
Adet = 50x50 pm . Assuming a blackbody source, St ean be taken to be unity. Noise 
equivalent power is evaluated by first eonsidering the various sources of noise within the 
deteetor. In microbolometers, the predominant noise source at lower readout bias 
eurrents is Johnson’s [55]. At room temperature, the RMS voltage of this noise 
eontribution is given by [57]; 

Vj^^{4k,TRAn^" ^7.0^V ( 22 ) 
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where R = 10 k^2 is the equivalent resistanee for noise analysis when the mierobolometer 
and 20-kf2 load resistor are serially eonneeted in the readout eireuit [55]. The frequeney 
bandwidth, Af = 0.29 MHz, is obtained from the 30 Hz frame rate at which the FPA’s 
160x120 pixel array is individually sampled. It follows that detector NEP, given by the 
ratio of VjN to detector responsivity, = 2x10^ V/W [55], is approximately 35 pW. 

As a result of optical modifications intended to maximize the amount of terahertz 
radiation received by the FPA, two of the optical parameters in Eqn. 4 (T/, Eno) differ 
from those for the original system. Early experiments using the camera’s stock 
germanium (Ge) lens indicated that the incident terahertz beam was strongly attenuated 
by anti-reflection coatings applied to the lens; to correct this deficiency, the Ge lens was 
replaced with a 1-inch diameter, 20-mm focal length bi-convex f/1 lens made of 
Tsurupica (PPE-l”-20mm-BC, Microtech Instruments). Eor focusing of the incident 
collimated beam, the lens was mounted 20 mm from the EPA, yielding an effective 
camera field of view of 53°. Tsurupica (formerly known by the trade name Picarin), 
measured to have a transmissivity (T/) of 0.65 at 2.8 THz, was also used as the source 
material for the 2-inch diameter, 4-mm thick cryostat window. The emissivity metric s, 
which is directly proportional to the membrane’s absorption of incident radiation, has not 
been well characterized for terahertz frequencies. By using infrared emissivity values, 
however, one can obtain an upper-limit of detector performance and establish whether 
external illumination is necessary under the most liberal of assumptions. Using s = 0.8 (a 
well-documented value for the 8-12 qm wavelength range [55]), the NETD of the 
mierobolometer is found to be approximately 3 K. This result, which is more than an 
order of magnitude higher than the commercially-specified NETD of the camera for 
passive infrared imaging (0.1 K), confirms the need for external illumination when using 
the mierobolometer camera at terahertz frequencies. 
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V. EXPERIMENTAL PROCEDURES 


A. IMAGING ARRANGEMENT 

Having completed the theoretieal and experimental foundations for development 
of a terahertz imaging system, next began the more vexing proeedure of aetually 
demonstrating that such an imaging technology is, in fact, achievable in praetiee. To this 
end, initial imaging trials were eondueted in 2005 by Karunasiri, Chamberlin and Robrish 
at Agilent Laboratories with the 3.6-THz QCL and IR-160 mierobolometer eamera as 
souree and sensor, respeetively [34], During these early feasibility trials, the QCL’s 
optieal eonfiguration was as shown in Fig. 13. To aoeommodate the stringent eooling 
eonditions required of the laser, the QCL was attaehed within a elosed-eyele refrigeration 
ehamber (maintained at 10 K temperature) using a eopper-based laser mount to aehieve 
good thermal eonduetanee. Also mounted to the eopper earrier were a one-ineh diameter, 
90-degree off-axis parabolie and an 80x80 mm flat mirror for eollimating and redireeting 
the laser beam through the eryostat window to the deteetion system. Both mirrors were 
gold-plated to allow maximum refleetion of the QCL beam; eopper was ehosen for the 
mounting material to optimize heat transfer away from the laser during operation. The 
mierobolometer was used in its original eommereial eonfiguration with the FPA sampled 
at a 30-Hz frame rate. 
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Fig. 13. QCL assembly used in initial imaging trials. QCL was mounted sueh that beam is 
emitted horizontally onto a flat mirror, then further refleeted off of a parabolie mirror before 
exiting the ehamber through a Tsurupiea-based window (from [37]). 


Initial imaging experiments conducted under this configuration failed to produce 
any signal whatsoever. Further examination suggested that this failure was likely due not 
to a lack of absorption by the pixel membrane layer, but, rather, the presence of 
antirefiective (AR) coatings on the camera’s original germanium-based lens element. 
Terahertz radiation is a full order of magnitude longer than the infrared radiation for 
which the microbolometer is designed; consequently, it followed that this AR coating 
was possibly acting as an absorbing medium to the QCL’s output beam. This suspicion 
was confirmed when, upon removing the lens element entirely from the IR-160, the 
terahertz beam was immediately detectable by the FPA as a (roughly circular) bright 
patch [34]. Camera output was recorded at a 18 FIz frame rate using frame grabber 
software (Fig. 14). 
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(C) (d) 


Fig. 14. Still images of video clips (taken by a frame grabber at 18 FIz frame rate) of real-time 
terahertz detection using a 3.6-TFIz QCL and IR-160 microbolometer FPA. (a) Unobstructed 
beam, (b) Beam with a metal pen passing between beam source and detector, (c) Beam 
obstructed by one sheet of bond paper, (d) Beam obstructed by two sheets of paper. 


Fig. 14(a) is a video recording of the imaging results produced using the virgin 
beam only; Fig. 14(b) is video of the same beam, disturbed by a metal pen passing 
through the region between the QCL and detector. Fig. 14(c) and Fig. 14(d), video 
recordings of the QCL beam obstructed by common bond paper, serve to illustrate the 
detection capability of the system under more demanding imaging conditions. These 
latter two images established the important fact that the QCL beam is still detectable after 
passage through one or more layers of opaque nonmetallic materials, thus confirming that 
this optical scheme could indeed be leveraged for use in practical imaging applications. 

Along with independent investigations conducted in the same year at the 
Massachusetts Institute of Technology (MIT) [35], these NPS-Agilent experiments were 
among the first to ever successfully demonstrate real-time, uncooled detection of 
terahertz radiation using a microbolometer FPA. Upon completion of initial studies. 
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Agilent researchers provided their 2.8 and 3.6 THz lasers, gratis, to the Naval 
Postgraduate School for purposes of continuing this research activity. At that time, 
Lowe, Karunasiri and the author began efforts to improve the detection results obtained 
using the same source/sensor system. Specifically, this entailed successfully evolving a 
system which merely detects terahertz radiation into one capable of high-resolution 
imaging of objects placed within the path of the beam. This clearly required the 
integration of a new lens element to replace the camera’s Ge-based lens, which had been 
found to frustrate initial attempts at terahertz detection. To this end, NPS purchased two 
Tsurupica-based lenses (one f/1 bi-convex lens of f=20 mm; one f/1.7 aspheric lens of 
f=50mm) from Microtech, Inc. that were specifically designed for use in optical 
systems employing terahertz radiation. Furthermore, an additional off-axis parabolic 
mirror was added to the optical path, external to the cold head, for steering and focusing 
the beam such that a collimated beam was presented to the region immediately in front of 
the camera. One of the two Tsurupica lenses was then placed one focal length away from 
the camera FPA to ensure a focused image on the focal plane. Finally, a Tsurupica-based 
window was installed in the cryostat to maximize transmission of the terahertz beam 
from the QCL through the cold head assembly. (Glass, which is known to partially 
attenuate terahertz radiation, was used for the original cryostat window.) Transmittance 
of Tsurupica in the terahertz regime was measured using Fourier Transform Infrared 
(FTIR) spectroscopy; as seen in Fig. 15, these measurements indicate that transmittance 
at 2.7 and 3.6 THz is approximately 65% and 61%, respectively. 
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Fig. 15. Measured transmittance of a Tsurupica window as a function of frequency in the 
terahertz regime. 


Imaging experiments were then attempted by plaeing various objeets (metallie 
and nonmetallie) between the external parabolie mirror and the IR-160 eamera. In spite 
of these optieal improvements, operation of the souree/sensor system with the two new 
lenses failed to produee any images that substantially improved upon the original results 
obtained at Agilent by Karunasiri, Chamberlin and Robrish. As before, metal objeets 
strongly attenuated the beam through absorption; nonmetallie objects were, in the main, 
effectively penetrated by the beam. However, attempts to focus the physical profiles of 
metallic and nonmetallie objects alike were decidedly unsuccessful [37]. 

Beginning in early 2007, the author and Karunasiri implemented further 
engineering changes in the optical configuration of the terahertz detection system. Two 
principal deficiencies were identified as possible reasons for the inability to produce 
improved imaging results. Firstly, two of the system’s four imaging elements (the one- 
inch diameter parabolic mirror and the 80x80 mm flat mirror) were internal to the 
cryostat and thus completely inaccessible (for adjustment) during imaging experiments. 
Secondly, the laser assembly was positioned approximately two inches from the cryostat 
window—a displacement which does not promote optimal transfer of laser energy from 
the cryostat. 
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To correct these engineering deficiencies, a new copper mount was designed and 
(with much-appreciated assistance from NPS/PH staff members George Jaksha and Sam 
Barone) fabrieated. As shown in Fig. 16, the new mounting bloek was mueh larger in 
size (to further improve the rate of heat transfer away from QCL during lasing) and 
contained no optical elements whatsoever. The QCL was mounted within a slot on the 
left side of the copper mount using low-temperature conducting epoxy. When enclosed 
within the eold head, this left side of the mount was separated from the cryostat wall by 
less than 1.5 em—thus allowing full passage of the QCL beam through the Tsurupica 
window. 



Fig. 16. Second QCL assembly design used within cryostat (with external cold head housing 
removed). Laser element is affixed to a copper carrier to facilitate transfer of thermal energy 
away from the laser. When lasing, beam is emitted outward, from QCL edge, along direction of 
red arrow (inset). 


The final optical arrangement was as shown in Fig. 17. External to the cryostat 
were plaeed a pair of gold-plated, 90-degree off-axis parabolic reflectors (fil and f/2, 
each 50.8 mm in diameter) for foeusing and steering the beam toward the eamera. To 
minimize lateral alignment error, these reflectors were mounted upon a single optical 
post, and were positioned such that the highly divergent beam exiting the cold head was 
focused to a eollimated state upon reaehing the microbolometer camera. To achieve this 
state, the f/2 GAP was displaced 25 mm from the f/1 GAP, which was positioned 32 mm 
from the eryostat window. The 20-mm focal length (f/1) Tsurupica lens was mounted to 


36 





the camera within a custom-designed lens assembly that replaced the original (discarded) 
germanium lens package mentioned previously. To properly focus the incident 
(collimated) beam, the lens was mounted one focal length (20 mm) from the FPA, 
yielding an effective camera field of view of 53°. Ideal illumination of the FPA was 
found to result for a camera position (as measured with respect to the centroid of the 
Tsurupica lens) with 90 mm displacement from the f/2 OAP. Imaging experiments were 
then performed by inserting various objects roughly midway between the two parabolic 
mirrors, which produced a focused image at the microbolometer FPA. As suggested by 
the spatially close integration of components in Fig. 17, the total optical path length of 
this new configuration was kept as short as possible (within focusing constraints) to 
minimize beam absorption by the ambient air. 



Fig. 17. Optical configuration used in most recent series of terahertz imaging experiments. 
Lower and upper mirrors (50.8 mm and 101.6 mm focal length, respectively) were used to focus 
and steer the terahertz beam emerging from the Tsumpica window of the cryostat to the focal 
plane array of the microbolometer (beam path is illustrated by red arrow). 


This new optical configuration immediately yielded a dramatic improvement in 

imaging results. These images, which are presented in detail in the following chapter, 

provided exceptionally high resolution of metal and non-metal materials alike under 

illumination by both the 2.7- and 3.6-THz QCL beams. 
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B. 


FTIR ANALYSIS 


Fourier transform infrared (FTIR) spectroscopy analysis (performed at a 
resolution of 1 cm'^ and a mirror speed of 0.64 m/s) provided important information 
about the spectral nature of the QCL beam. Perhaps most notably, it revealed that the 

“3.4-THz” QCL provided by Agilent was, in fact, a 3.6-THz laser. Apart from this 
minor but important technical revelation, FTIR spectroscopy established the spectral 
lineshape associated with a typical QCL beam. As seen in Fig. 18, the full-width, half¬ 
maximum (FWHM) spectral width of the QCL beam is approximately 45 GHz, with 
resonance occurring at 3.57 THz (with both values exhibiting slight sensitivity to applied 
duty cycle). As expected, QCL output power was found to generally increase with 
increasing duty cycle. FTIR analysis also provided an opportunity to further understand 
the relationship between duty cycle and output power over a range of applied bias 
currents. 



Frequency (THz) 


Fig. 18. Optical configuration used to conduct FTIR measurements of 3.6-THz QCL beam. A 
gold-coated, f/1 parabolic reflector was used with a silver-coated parabolic reflector for focusing 
and directing the beam into the external port of the FTIR spectrometer. Measurements were 
taken with the QCL operating under an applied bias of 1.6-1.7 A. 


It is perhaps worth clarifying, at this point, that this disparity (3.4-THz vs. 3.6-THz) is the reason 
that two conference presentations by the author in early 2007 contain “3.4-THz” in the title. These 
presentations were prepared under the assumption that, per manufacturer’s claims, 3.4-THz was indeed the 
lasing frequency of the QCL. Not until after the conclusion of the two conferences was the above- 
mentioned FTIR analysis performed, yielding the 3.6-THz value which is now known to be the actual 
resonance frequency. 
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Fig. 19 shows measured output power of the laser for applied bias of 1.0-2.2 A, as 
measured by FTIR spectroscopy. The laser was operated at duty cycles of 5-25% 
(170-830 ns pulse length) to determine optimal pulse conditions for use in imaging 
experiments. Signal strength was measured by feeding the QCL beam through the 
external port of the spectrometer and measuring the height of the resulting peak. Higher 
duty cycles were found to generally produce higher output power; however, as further 
investigations confirmed, such conditions tend to rapidly warm the laser assembly via 
Joule heating effects. Each FTIR measurement (consisting of 16 scans) was taken over 
approximately 35 seconds with an initial cryostat temperature of lOK; therefore, the 
results presented in Fig. 19 are only representative of the initial output power from the 
laser. 



Fig. 19. Output power as a function of current for various duty cycles (D). Each data point 
was collected by FTIR spectroscopy using 16 scans at a resolution of 1 cm'* and a mirror speed of 
0.64 m/s. 
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c. 


FPA SENSITIVITY 


Another, equally important method of QCL beam eharaeterization involves 
assessing FPA sensitivity to the ineident beam. To the extent that the mierobolometer 
eannot produee images with a beam that is poorly absorbed by the FPA, this evaluation 
was essential to optimizing the image quality obtainable by the souree/sensor system. 
This was done by measuring FPA signal intensity, for a speeifie duty eyele, as a funetion 
of the laser bias eurrent. Fig. 20 reports the measured signal intensity of the laser while 
operating at a 15% duty eyele, overlaid with a eorresponding V-I eurve. Signal intensity 
of eaeh mierobolometer image was ealeulated by using a eomputer algorithm to 
deeonstruet the image, evaluate the brightness of eaeh of the 19,200 pixels, and sum the 
result. For these partieular operating eonditions, threshold lasing was aehieved at 
approximately 1.0 A; lasing was found to eease beyond 2.2 A primarily due to (1) 
misalignment of energy levels due to large bias in the deviee and (2) exeessive heating of 
the aetive region, originating from this same bias eurrent, and resulting in thermal 
population of the middle laser level. Maximum beam intensity was aehieved at 1.8 A 
with a bias voltage of about 10.5 V. For illustrative purposes, representative images of 
FPA output are ineluded as inset frames within Fig. 20. The dual nature of this figure 
(reporting both voltage and intensity vs. bias eurrent) is useful in that it illustrates a 
reeurring trend that is observable for this and other duty eyeles. Speeilieally, measured 
intensity tends to rapidly deeline for bias eurrents above that value at whieh the QCL V-I 
eurve eeases to be linear (approximately 1.9 A). 
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Fig. 20. Measured voltage aeross the laser and normalized image intensity of the QCL, 
operated at 15% duty eyele, as a function of current. The image quality associated with specific 
levels of signal intensity is illustrated by the inset images, taken by microbolometer FPA, of an 
unobstructed beam at operating currents of 1.0 A, 1.3 A, 1.8 A, 2.1 A, and 2.2 A. 


D. IMPACT OF QCL TEMPERATURE ON IMAGING RESULTS 

As alluded to above, the preceding graphs representing output power and signal 
intensity must be viewed with the caveat that they are, in the main, only indicative of 
beam strength over relatively short periods of time. Especially at higher bias currents 
and/or higher duty cycles, they reflect values obtained before Joule heating effects 
accumulate to the point of adversely affecting QCL performance. These effects are 
certainly not negligible; rather, they figure prominently into the system’s practical 
viability as an imaging apparatus. Fig. 21 serves to illustrate cryostat temperature, at 
various duty cycles, over prolonged periods of time. Notwithstanding the fact that the 
cryostat was set to 10 K in each case, close-cycle refrigeration is not sufficiently robust to 
completely overcome Joule heating of the laser during operation. As a result, 
temperature of the laser assembly rises until reaching a steady-state value that is 
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dependent upon the applied voltage bias and selected duty cycle. (For the conditions 
under which Fig. 21 data were collected, these temperatures were 18, 21, 24, and 32 K 
for duty cycles of 10%, 15%, 20%, and 25%, respectively). The only case in which 
temperature equilibrates at a value appreciably below its peak value is for the lowest duty 
cycle under examination (D = 10%). Herein can be clearly seen the competing effects of 
Joule heating and active cooling by the cold-cycle refrigerator. For this particular case 
(in which heating effects are more modest), the initial rapid temperature rise is eventually 
mitigated by the much slower cooling process. This important observation suggests that, 
for low duty cycles, FPA-detected signal strength may actually improve with increasing 
time after an initial reduction that is transient in nature (albeit not to the high signal level 
obtained when lasing is first initiated). At the other extreme (D = 25%), cryostat 
temperature was found to monotonically increase for an appreciably longer period of time 
than lower duty cycles before reaching a steady-state value. Indeed, the time required to 
achieve thermal equilibrium was found to rise proportionally with the duty cycle under 
which the laser was operated. 



Duration of Lasing (min) 

Fig. 21. Effects of duty cycle on laser performance: Signal intensity detected by FPA while 
operating QCL, at 300 kHz pulse rate and applied bias of 9.7-10.5 V, over extended period of 
time. (Cryostat set to 10 K.) 
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As predicted, these observed temperature changes are, in turn, responsible for a 
corresponding decline in QCL signal strength. Again, this is likely due to the reduced 
population inversion that results from thermal excitation of electrons between the 
quantized energy states of the laser. Fig. 22 quantifies the extent to which duty cycle 
affects laser beam intensity over extended periods of time. From these results, it is clear 
that QCL operation at a 25% duty cycle is fundamentally different from that produced 
using lower duty cycles. Temperature rises occur under each scenario. But, whereas 
signal intensity remains roughly constant in time for duty cycles of 10-20%, it drops 
precipitously (and immediately) for values above this range. In light of these 
observations, it can be posited that there is a threshold value of duty cycle, between 20% 
and 25%, above which the laser is no longer capable of serving as a useful terahertz 
imaging illumination source for more than a brief period of time. 



Fig. 22. Effects of temperature on laser performance for various duty cycles (D): Signal 
intensity detected by FPA while operating QCL, at 300 kFlz pulse rate and applied bias of 9.7- 
10.5 V, over extended period of time. (Cryostat set to 10 K.) 
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These fundamentally distinet results ean be explained by eonsidering the 
magnitude of temperature ehange elieited by eaeh seenario. At lower duty eyeles, the 
relationship between eryostat temperature and duty eycle was found to be approximately 
linear (with slope of 0.6K/peroent). For duty eyeles above 20%, though, linearity of the 
relationship eollapses and temperature of the laser assembly rises dramatieally—nearly 
extinguishing the output beam entirely. As mentioned previously, this is likely due to 
thermally-aetivated eleetrons populating the lower laser level from the first eonduetion 
subband. During lasing at low temperatures, the lower laser level is depopulated via 
photon-phonon interaetions; however, at elevated temperatures, this proeess is 
frustrated—^preventing population inversion between the upper and lower laser levels. 

In light of these findings, imaging experiments were eondueted under different 
operating eonditions, depending upon the type of deteetion desired. For single-frame 
images taken immediately after the eommeneement of lasing, large voltage bias and high 
duty eyele were used to produee a high-power (albeit short-lived) output beam. In 
general, though, these aggressive eonditions were not found to produee an exeeptionally 
superior output power. Furthermore, they would have little utility in any praetieal 
imaging system. As a result, most imaging—for both still frames and extended video 
reeording—^was performed under more modest operating eonditions (typieally 10-11 V 
applied bias, with 15-20% duty eyele.) 

Fig. 23 provides further insight into the temporal behavior of QCL temperature 
during typieal operation. When the QCL is first turned on with the applieation of an 
eleotrieal bias (1.7 A in the present ease), temperature rises, in logarithmie fashion, to a 
steady-state value over the eourse of several minutes. After removing the applied voltage 
bias, the sensor temperature drops rapidly (initially in linear fashion, at a rate of 
approximately 0.2 K/see). Within roughly three minutes, the deviee returns to the 10 K 
temperature at which the chamber thermostat is set. 
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Fig. 23. Cryostat temperature changes arising from operation of the QCL at 20% duty cycle 
and applied bias of 10.4 V (1.7 A). Under these conditions, the cryostat reaches a steady-state 
temperature of 29 K after approximately eight minutes. Upon ceasing current flow to the laser, 
cryostat temperature drops at an approximate rate of 0.2 K/sec. 
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VI. IMAGING RESULTS 


The two QCLs used in these imaging experiments follow the same basic design 
principles and were fabricated using similar methods; however, the quality and 
persistence of imaging results yielded by the lasers differ dramatically. In general, results 
obtained using the 2.7-THz QCL as an illumination source were inferior in all respects to 
those produced using the 3.6-THz laser. Much of this deficiency appears to be directly 
attributable to the effects of elevated temperature on the laser, even though the rate of 
temperature rise for the two devices were roughly the same (from lOK to ~30K in 3-5 
minutes). For instance, under nominal operating conditions, the 2.7-THz laser produced 
high-quality video imaging for a duration of 4-6 seconds; longer intervals of time rapidly 
produced a “wash-out” effect in which the beam was almost entirely extinguished and 
undetectable by the microbolometer system. The beam was again reproducible by cutting 
power from the device and waiting 5-10 minutes for the Joule heating in the laser to 
dissipate. Conversely, the 3.6-THz laser was capable of producing a consistent, high- 
intensity beam for extended periods of time (10 minutes or more). 

Much of this disparity in the performance of the two laser devices can be 
understood by considering the relative energy gap produced by the multiple quantum well 
structure in the laser active region. As previously explained, both lasers operate off of a 
three-level lasing scheme within the conduction band; the effective energy gap exploited 
by the 2.7 THz laser and the 3.6 THz laser are 11 meV and 15 meV, respectively. This 
energy gap difference, which is more than 35% of the 11 meV value, is particularly 
significant when considering the fact that thermal population of the upper lasing level is 
exponentially increasing with temperature. At 300K, the Boltzmann factor for the 
2.7-THz and 3.6-THz lasers are roughly equal, giving a ratio of electron populations for 
levels 3 and 2 (N3/N2) of about 0.65 and 0.56, respectively. At cryogenic temperatures, 
however, the difference in Boltzmann factor becomes much larger. At 10 K, for instance, 
N3/N2 for the 2.7-THz laser is more than a factor of 75 greater than that of the 3.6-THz 
QCL. 
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Imaging data produced using the 2.7-THz and 3.6-THz lasers demonstrate the 
signifieant value of having more than one terahertz laser available for this series of 
experiments. These data suggest that QCLs designed with higher resonant frequencies 
ean, in general, be expeeted to produee superior illumination eonditions for terahertz 
imaging. As corroborated by the quality of still and video images in the following two 
seetions, this improvement is manifest in two ways. The most notable and important 
improvement is in the persistence of illumination during extended video imaging. 
However, it is evident that the higher-frequeney QCL provides a slightly superior quality 
of single-frame images as well. 

A. IMAGING AT 2.7 THz 

Using the optieal arrangement described above, imaging experiments were 
performed by inserting metallic objects, wrapped within various obseurants (plastic, 
paper, and cloth), roughly midway between the two parabolie mirrors sueh that a focused 
image fell on the mierobolometer FPA. For single-frame imaging, the laser was typically 
operated at a 300 kHz repetition rate, a bias of 1.3-1.6 A, and 15% duty eyele. Fig. 24 is 
representative of the results obtained and demonstrates the methodology used to refine 
the image quality of a utility knife blade wrapped in opaque plastie tape (Fig. 24a). (One 
edge of obscuring tape was eut with a non-linear contour and notched so as to better 
identify orientation of the knife and facilitate evaluation of the system’s imaging 
success.) 
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(c) (d) 


Fig. 24. Imaging of a small utility knife blade wrapped in opaque plastic tape, (a) 
Conventional digital photograph. Red dotted region represents approximate area of illumination, 
(b) Single frame image of blade assembly illuminated with 2.7-THz QCL radiation and imaged 
with microbolometer camera, (c) Image generated by computationally averaging 50 individual 
frames, (d) Fifty-frame composite image, refined using MATLAB image processing utility 
software. Orphan particle is visible within red circle. 


Although single-frame images showed good contrast between the two objects 
(Fig. 24b), further improvements were achieved by taking multiple individual images of 
the stationary ensemble. When image averaging techniques are used to combine these 
frames into a composite image (Fig. 24c), noise effects are averaged out and overall 
clarity is improved. Finally, using noise reduction software applications, still greater 
contrast can be achieved (Fig. 24d). Fig. 24d was obtained by operating upon Fig. 24c 
with a series of MATLAB image processing algorithms; contrast-limited adaptive 
histogram equalization (CLAHE), decorrelation stretch with 5% tolerance saturation, and 
intensity adjustment with 1% saturation. The same imaging methodology is used for 
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imaging of a large paperclip in Fig. 25. Fig. 26 shows imaging results for three 
additional objects under similar optical conditions. 



(c) (d) 


Fig. 25. Imaging of a large paper elip wrapped in opaque plastic tape, (a) Conventional 
digital photograph. Red dotted region represents approximate area of illumination, (b) Single 
frame image of blade assembly illuminated with 2.7-THz QCL radiation and imaged with 
microbolometer camera, (c) Image generated by computationally averaging 50 individual frames, 
(d) Fifty-frame composite image, refined using MATLAB image processing utility software. 
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(d) (c) (0 


Fig. 26. Various objects, wrapped within two layers of opaque plastic tape, imaged using 
2.7-TFIz QCL and microbolometer FPA. Upper images are conventional photographs taken with 
a compact digital camera; lower images are 50-ffame averages of terahertz images enhanced with 
MATLAB noise reduction post-processing, (a, d). Dentist’s pick, (b, e). Dissection scalpel with 
curved blade, (c, f). Plastic tie bent into a loop. 


The advantage offered by image averaging and/or eomputer-based post- 
proeessing techniques is perhaps best appreciated by noting the small speck of 
(nonmetallic) foreign material which, having become accidentally embedded within the 
plastic tape, can be seen in terahertz images near the top of the illumination region 
(circled in Fig. 24d). In the unmodified single-frame terahertz image (Fig. 24b), the 
contaminant is barely distinguishable within the image’s background noise. In Fig. 24c 
and Fig. 24d, however, it is clearly recognizable as a distinct object. 

As previously mentioned, the primary limiting factor to this imaging arrangement 
is rapid heating of the laser due to Joule heating from the input current. Heating effects 
were mitigated by chopping the 300-kHz driving signal with a 20-Hz gating frequency 
and slow sweep rate to allow the device to cool intermittently. Due to its intermittent 
nature, such an illumination scheme is not useful for extended-duration imaging of 
moving objects. To allow imaging for several continuous seconds, the external gate 
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generator was disabled and input eurrent and duty eycle of the input signal were both 
appreeiably deereased (to roughly 0.93 A and 8%, respeetively) to mitigate heating 
effects. 

The results of video imaging with the 2.7-THz QCL are presented, in .avi format, 
in Fig. 27. Despite the resulting reduction in average laser power, real-time video recordings 
under these operating parameters are only slightly inferior to still-frame images. Although raw 
video was of acceptable quality, video filters (blur and box blur; sharpen, level and HSV 
adjustments) were applied to each of the four images to improve overall contrast. 



(c) (d) 


Fig. 27. Still images of video clips (with post-processing filters applied) from imaging 
experiments conducted using microbolometer FPA and 2.7-TFlz QCL operated at a 300 kFlz 
pulse repetition rate, 8% duty cycle, and applied bias of 1.0 A. All objects were wrapped in two 
layers of opaque plastic tape, (a) Steel utility blade, (b) Dentist’s pick, (c) Large paperclip, (d) 
Plastic tie bent into the shape of a wire loop. 


Significant diffraction effects, produced by the influence of various system optics on the 
terahertz beam (QCL aperture diffraction, as well as etalon effects arising within the 
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Tsurupica lens and cryostat window), are apparent in the form of sharply-defined, 
concentric fringe patterns appearing throughout each of the terahertz images. 


B. IMAGING AT 3.6 THz 

Like those conducted at 2.7-THz, imaging experiments at 3.6-THz were 
performed with a 300-kHz PRF and a 30-Hz FPA frame rate. Due to its higher resonant 
frequency (and the correspondingly larger energy gap used for stimulated emission), 
however, the 3.6-THz QCL was able to tolerate higher temperatures and therefore operate 
effectively under more robust electrical conditions. For still-frame and extended video 
imaging alike, the laser was nominally operated with an applied bias of 1.6-1.9 A and a 
15-20% duty cycle. 

Fig. 28 is a collection of video clips demonstrating the type of real-time imaging 
that is obtainable under the aforementioned optical arrangement. All images were taken 
(using commercial frame grabber software) by the IR-160 microbolometer camera, at a 
30-Hz frame rate, under illumination by the 3.6-THz QCL. The QCL was operated with 
a bias of 1.9 A at a 300-kHz pulse rate and 20% duty cycle. 
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(c) (d) 


Fig. 28. Still images of video clips (with post-processing filters applied) from imaging 
experiments conducted using microbolometer FPA and 3.6-TFIz QCL operated at a 300 kFIz 
pulse repetition rate, 20% duty cycle, and applied bias of 1.9 A. (a) Steel utility blade obscured 
by two layers of opaque plastic tape, (b) Utility blade concealed with common bond paper, (c) 
Mylar® film cut in the shape of the Ironman Triathlon® “M-Dot” logo and enclosed by two layers 
of opaque plastic tape, (d) Finger of a black polyurethane glove containing an Allen wrench, 
which in turn is wrapped with a small swatch of opaque plastic tape. 


Interest in terahertz imaging is largely driven by the fact that, whereas terahertz 
radiation is strongly reflected by metals, it penetrates most non-metallic materials with 
little-to-no attenuation. This provides good image contrast between metals and 
obscurants of non-metallic constitution. As is seen in Fig. 28, the contrast is particularly 
striking in the case of plastic materials. Steel blades, such as the one imaged in 
Fig. 28(a), can be clearly seen when wrapped in opaque plastic tape—and remain visible 
even when enclosed by as many as eight layers of tape. To a lesser extent, metals are 
also detectable through cloth and common bond paper (Fig. 28(b)). Biaxially-oriented 
polyethylene terephthalate (boPET) film (more commonly known by the trade name 
Mylar®) was found to also strongly absorb terahertz radiation—even when used in 
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exceptionally thin quantities (Fig. 28(c)). Finally, it is noteworthy that this terahertz 
imaging scheme is also capable of imaging non-metallic materials within other non- 
metallic materials. In Fig. 28(d), the terahertz beam is sufficiently attenuated by both 
polyethylene and plastic tape that each material can be clearly distinguished from the 
other. 

Fig. 29 contains another set of four images, taken under the same illumination 
scheme as in Fig. 28. Herein, more challenging sets of imaging conditions were 
investigated—demonstrating that available power from the QCL is somewhat inadequate 
for penetrating thicker obscurants. However, image resolution is still sufficiently robust 
so as to allow identification of the various objects. Furthermore, it is in such demanding 
imaging arrangements that the utility of image processing and/or color addition becomes 
apparent.tt In each of the four imaging samples represented within Fig. 29, video flies 
were decomposed using VirtualDub (a video capture/processing application available for 
free use under the GNU General Public License) into a series of individual images; these 
images were then adjusted for contrast, brightness, and saturation level using a batch 
operation in Microsoft® Picture Manager to obtain optimal differentiation of the distinct 
objects under illumination. Final reassembly of the images was performed using 
VirtualDub. By employing this approach (video decomposition, image adjustment, and 
reconstitution), one obtains an improved level of editorial control over the quality and 
resolution of the original video file. Nevertheless, it should be emphasized that the final 
imaging product does not represent a dramatic improvement over the original video 
results; direct FPA video output still allows full differentiation of the objects of interest. 

Fig. 29(a) is the now-familiar utility blade from previous images concealed 
behind a single layer of relatively thick (2 mm) Plexiglas. Fig. 29(b) demonstrates the 
strong terahertz attenuation associated with even thin layers of sponge-like foam packing 
material. In this image, an embedded 0.5-mm lead from a mechanical pencil is barely 
resolvable. Fig. 29(c) is extended video recording of terahertz imaging of a U.S. $5 bill. 


tt For the remainder of this thesis, for purposes of easy identification and categorization, red 
colorization is used in images of objects concealed within obscurants; blue colorization is used in images of 
paper-based materials directly illuminated by TFlz radiation. 
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Notable in this image are two types of government seeurity (anti-eounterfeit) features: a 
vertieally-oriented narrow plastie thread (~ 1 mm in width) and an internal watermark of 
President Abraham Lincoln (3 mm long by 2 mm wide) on the left and right extremes of 
the bill, respectively. (Neither security feature is visible under normal observation; 
however, each can be detected when back-illuminated by a white light source.) 



(c) (d) 


Fig. 29. Still images of video clips (with colorizing and post-processing filters applied) from 
imaging experiments conducted using microbolometer FPA and 3.6-TFlz QCL operated at a 300 
kFlz pulse repetition rate, 20% duty cycle, and applied bias of 1.9 A. (a) Steel utility blade 
obscured by a layer of Plexiglas, (b) mechanical pencil lead embedded in foam, (c) Watermark 
of President Abraham Lincoln from a U.S. $5 bill (embedded video clip also includes imaging of 
plastic security thread), (d) “NPS” (“N” in still-frame image) written in number-two pencil, on a 
sheet of common bond paper. 


Great Britain Pound (GBP) notes arguably contain some of the most robust 
security features of any currency in the world. For this (and other) reasons, it is also a 
superb choice for examination under terahertz imaging. For purposes of this research, a 
GPB £5 note was used as the imaging specimen. As with U.S. currency notes, British 
notes contain a watermark which (when using a conventional white light source) is only 
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visible under back-illumination (Fig. 30a). As seen in Fig. 30b, this watermark is clearly 
discemable in terahertz images as well. While the terahertz image does not contain the 
same level of detail and clarity produced by rear-illumination of visible light, it does 
potentially provide an additional layer of security against counterfeit practices by making 
watermark duplication more difficult (insofar as any attempt at watermark forgery would 
need to be visible under both imaging schemes). 



(c) (d) 


Fig. 30. £5 Great Britain Pound (GBP) note, (a) Watermark image of Queen Elizabeth II, as 
seen under illumination by forward-incident and rear-incident white light source (top-left and 
bottom-right comers, respectively). Only under the latter scheme is the watermark visible, (b) 
Terahertz imaging of watermark using microbolometer FPA and 3.6-TFIz QCL operated at a 
300 kHz pulse repetition rate, 20% duty cycle, and applied bias of 1.9 A (with post-processing 
filters applied), (c) £5 security thread as seen under illumination by forward-incident and rear- 
incident white light source (top-left and bottom-right corners, respectively), (d) Terahertz 
imaging of security thread (using same operating parameters as in (b)). Notable in this animation 
is the elaborate cross-thread stmcture, a security feature which is relatively undetectable using 
white light sources. 


More impressive is the ability of terahertz imaging to detect security features 

within GBP notes which are relatively undetectable under conventional white light 

illumination. British notes contain an internal security thread which is decidedly more 
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complex than U.S. currency. The thread (Fig. 30b), which appears as a dashed and solid 
line when viewed under illumination by forward-incident and rear-incident white light 
sources, respectively, also contains an elaborate cross-thread structure which is not 
immediately apparent, even with the highest levels of scrutiny, under visible light. As 
seen in Fig. 30d, terahertz imaging reveals that these cross-threads, which are remarkably 
elongated in the lateral dimension, appear at regular intervals along the entire length of 
the security thread. (The strong attenuation of terahertz radiation by these cross-stitches 
suggests that they are quite likely metallic in nature.) Were terahertz imaging to be 
adopted as a countermeasure against GBP forgery, an aspiring counterfeiter would be 
forced to produce internal threads which contain the same unique optical properties as 
this material (solid straight line under visible light; stitching under terahertz radiation)—a 
difficult prospect indeed. It is worth noting that the existence of this cross-thread feature 
is, to the best of the author’s knowledge, relatively obscure. No mention of it can be 
found in the public domain, even within banking and government web sites which 
purport to document the complete set of security features (some of which are only 
viewable under black or ultraviolet light) that are embedded into each nation’s currencies. 

C. IMAGING AT 0.7-3.0 THz USING INTRACAVITY DFG 

In the very final stages of this dissertation research (immediately after concluding 
a presentation at the 2008 Photonics West conference), the author was approached by 
Stanford University research professor Dr. Konstantin Vodopyanov with a very timely 
and promising opportunity for collaboration. Dr. Vodopyanov’s research group had 
recently developed a specific lasing scheme in which milliwatt-scale terahertz radiation is 
produced by applying difference frequency generation (DFG) techniques to a 1064-nm 
Nd;YV 04 source laser [58]. Through the introduction of a quasi-phase-matched 
nonlinear crystal into the beam path (for the particular research effort cited here, 
periodically-poled lithium niobate (LiNb 03 ) or GaAs crystals are used), an optical 
parametric oscillator (OPO) mechanism is produced in which the Nd:YV 04 pump beam 
(with radial frequency of C 0 p=l.77 x lO'^ rads/sec) can be down-converted into two 
approximately-identical beams, of radial frequency cOg and coi (corresponding to the signal 
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and idler beam, respeetively). Plainly, eonservation of energy requires that Wi + ©i = ©p. 
By judieious seleetion of the signal and idler frequeneies, then, any frequeney below that 
of the pump laser ean (in theory) be produeed by DFG of the two slave beams using 
nonlinear erystals. For the generation of terahertz radiation, the only eritieal requirement 
is that the signal and idler beams differ in frequeney by that value desired for the DFG- 
produeed beam. The speeifie proeess eurrently in use by Dr. Vodopyanov’s group is 
doubly-resonant (both the signal and idle beams are produeed and amplified within the 
eavity), near-degenerate (©s and ©i are almost, but not quite, identieal) parametrie 
amplifieation. The partieular frequeney produeed by the DFG seheme is determined 
solely by the refraetive index of the nonlinear erystal, whieh ean be eontrolled through 
adjustments in erystal temperature (these temperature shifts are produeed through 
eonduetion effeets arising from heating of the erystal on a thermal stage). 

The proposed eollaboration was realized on 29 Feb 2008, when NPS and Stanford 
researehers met and worked together in Ginzton Laboratory to eonduet experiments 
involving the OPO and the IR-160 mierobolometer eamera. In these imaging trials, the 
IR-160 eamera replaeed the deuterated triglyeine sulfate (DTGS) deteetor in the optieal 
arrangement and used for deteetion of 3.02, 2.8, 2, and 0.7 THz output beams. Due to the 
frequeney-mixing seheme employed, there is a strong infrared eomponent (attributable to 
the ~2 pm signal and idler beams) in the laser’s virgin output signal. This eomponent 
was effeetively eliminated by ineluding a Mieroteeh low-pass filter, whieh has -100% 
attenuation for the infrared eomponents involved, into the beam path. The IR-160 was 
mounted (without the 20-mm Tsurupiea lens used in the previous experiments) one foeal 
length from a 50-mm Tsurupiea lens used in the Stanford optieal system. As seen in Fig. 
31 and Fig. 32, these beams were imaged sueeessfully. Although Fig. 31 is dominated by 
infrared eontributions from the signal and idler beams, it does provide valuable insight 
into the quality of beam used for generation of the terahertz signal. 
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Fig. 31. Three-dimensional intensity map of the Stanford OPO laser beam during production 
of 2.8 TFIz radiation, as measured by IR-160 focal plane array (without low-pass filter). Due to 
the frequency-mixing scheme, beam intensity is predominantly attributable to the two ~2 pm 
signal and idler beams. 



Fig. 32. Colorized frame capture of IR-160 camera under illumination by 2.8 THz OPO 
beam (with Microtech lowpass terahertz filter included in beam path.) 


The most noteworthy—if not surprising—produet of these experiments involves 
examination of the 0.7-THz beam. Despite a relatively low signal-to-noise ratio, the 
0.7 THz beam was sueeessfully deteeted and imaged by the IR-160 eamera. To the 
author’s knowledge, this is the lowest frequeney of radiation ever sueeessfully deteeted 
using mierobolometer technology—and suggests that microbolometer pixel membranes 
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remain absorptive at frequencies much lower than those for which they were designed. 
While confidence is high that the detected signal is indeed (and exclusively) the 0.7-THz 
radiation (and not the result of infrared beams bleeding through the low-pass filter), 
follow-up experiments are being designed that will enable the researchers to definitively 
confirm this conclusion. 
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VII. POTENTIAL IMPROVEMENTS TO THE IMAGING SYSTEM 


Having established the fundamental feasibility of real-time terahertz imaging 
using satisfaetory but far-from-ideal eomponents, the Sensors Researeh Group (SRG) has 
reeently made sizable inroads towards further improving both the quality and 
eompaetness of a terahertz imaging system. This task has been approaehed from both 
souree and the sensor perspeetives. These inelude, most prominently, pursuing 
proeurement of another, higher-power QCL from Jerome Faist’s researeh group at the 
University of Neuehatel, as well as the previously-mentioned eollaborative opportunity 
with Stanford University. While average output power levels of the Stanford laser are 
eurrently inferior to the NFS QCL system, the Stanford laser exhibits remarkable 
tunability over terahertz frequeneies (0.7-4.8 THz), and produees an output beam that is 
impressively “elean” (spot size is almost perfeetly eireular, and is not blemished by the 
fringing effeets that are so pernieious in the NFS QCL). This latter eharaeteristie is likely 
to be partieularly valuable for imaging applieations—espeeially for objeets exhibiting 
fine detail, sueh as watermarks embedded within United States dollar (USD) or GBF 
notes. Fringing effeets are little more than a minor liability when imaging bulk objeets; 
however, images of objeets eontaining print or fine detail are severely degraded by these 
diffraetion-limited effeets. Furthermore, the broad tunability of the Stanford laser 
provides an opportunity to more elosely examine threshold deteetion limits of the 
mierobolometer eamera (determining high- and low-frequeney eutoffs for the eamera’s 
sensitivity). Finally, the Stanford laser operates sueeessfully under room temperature 
eonditions, obviating the eurrent need for bulky and time-eonsuming eryogenic eooling 
equipment/proeesses and furthering the goal of developing a eomplete system (as 
opposed to only the sensor) that is eapable of operating in a truly uneooled state. 

A. DIFFRACTION EFFECTS 

While not critically debilitating, the manifestation of diffraetion effeets acts to 
degrade overall resolution of the terahertz imaging system—particularly in the case of 
objects containing fine detail (e.g., images and/or text). Indeed, these effects currently 
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serve to elassify the imaging seheme as a diffraetion-limited system. These diffraetion 
effeets are prineipally driven by (1) the long wavelengths assoeiated with terahertz 
radiation (-100 pm), (2) the extraordinarily small output aperture of the QCL 
(14 pm X 200 pm), and the presenee of multiple optieal eomponents along the full extent 
of the beam path. As the first of these eonditions is, by definition, an unavoidable 
eonsequenee of terahertz radiation, attention must be direeted to eonsideration of the 
latter two eonditions. Small aperture size is an inherent feature of QCLs and thus not 
eontrollable under the eurrent imaging seheme; however, adoption of a new laser souree 
that operates upon entirely distinet physieal laws is one straightforward and reasonably 
aeeessible solution to this teehnieal issue (a point whieh argues strongly in favor of 
aggressively pursuing eollaboration opportunities with Stanford). The last of these 
eonditions, the influenee of opties elements, is eurrently being remedied through the 
proeurement of a wedged eryostat window. By replaeing the uniformly-thiek (4-mm) 
window with a wedged window of 3° inelination, any fringes whieh arise due to etalon 
effeets within the flat window should be eompletely eliminated. Eliminating diffraetion 
effeets assoeiated with the other three lensing media is a more intraetable prospeet; 
experimentation has demonstrated that they are, in large measure, absolutely neeessary 
for produeing an aeeeptably foeused image at the FPA. However, to the extent that the 
Tsurupiea lens has been found to be responsible for a large proportion of the eumulative 
fringing phenomena, it is worthwhile to examine Tsurupiea lenses of other diameters and 
foeal lengths. These options are eurrently being investigated. 

B. OUTPUT POWER 

As is often the ease with basie researeh, the prineipal goal of this dissertation 
work is one of feasibility demonstration; a theoretieal and experimental demonstration 
that real-time terahertz imaging ean be sueeessfully performed with the use of uneooled 
mierobolometer teehnology. Within this set of expeetations, the results obtained here are 
both sueeessful and promising—and offer solid proof that sueh an imaging approaeh is 
worthy of further investigation. However, it is undeniable that the low output laser power 
from the eurrent terahertz souree plaees rather stringent eonstraints upon the thiekness 
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and nature of materials through whieh imaging ean be performed—and eonstrains the 
eonditions under whieh the mierobolometer’s imaging eapabilities ean be tested. With 
less than one milliwatt of average power at 3.6 THz, metal eoneealed by approximately 
4-mm of plastie opaque tape ean be sueeessfully imaged. But terahertz imaging will not 
be eompetitive with other seeurity imaging teehniques (e.g., x-ray baekseatter) until it is 
eapable of diseeming objeets through and within a wide range of thiek non-metallie 
materials. Terahertz transpareney for many sueh materials has already been 
demonstrated; however, even relatively minor eleetromagnetie attenuation ean and does 
serve to dramatieally deerease the amount of terahertz power falling upon the FPA. A 
higher laser output power is eertain to inerease the thiekness (and quite likely the breadth) 
of materials through whieh metallie objects can be imaged. The SRG is currently 
exploring (through correspondence with Prof Jerome Faist, the acknowledged pioneer 
(and co-inventor) of QCL technology) the prospect of obtaining another, higher-power 
QCL for use in future research experiments. 

Another potential solution to the problem of low laser power lies with 
transitioning to an entirely different laser device. There exists an important experimental 
niche that can be filled by a custom-designed, “tabletop”-size free electron laser (FEL). 
Despite the fact that various U.S. laboratories have free electron laser facilities that are 
available for use by academia, most such lasers do not have the pulse characteristics that 
allow imaging using microbolometer camera systems. Many systems are not capable of 
lasing in the terahertz regime; others are capable of high-power terahertz pulses, but with 
slow repetition rates and exceptionally short pulse widths that are not amenable to 
microbolometer detection. 

An example serves to illustrate the problem associated with these pulse 
characteristics. The far-infrared FEL system at El.C. Santa Barbara is capable of lasing at 
63-338 pm with a peak power of 1-6 kW per pulse [59]. With a specified pulse length of 
1-20 ps and a pulse repetition frequency (PRE) of 0-7.5 Hz, the absolute upper limit of 
laser average output power is given by: 

Pa., = P,ea,-t.uise-fpRF = (6 ^1F)(20/is)(7.5) = OOOmlF (23) 
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Although this average output power is roughly two orders of magnitude greater 
than that eurrently produced by the aforementioned QCLs, the characteristics of the pulse 
itself are almost entirely incompatible with real-time microbolometer detection. The 
device’s voltage source, which is essentially an enormous Van de Graff generator, 
requires time to build a sufficient electric charge for each pulse; this is the limiting factor 
which is responsible for the slow PRF. Therefore, using this illumination source roughly 
reduces the effective frame rate of the camera fourfold: from 30 Hz to 7.5 Hz. 

More importantly, microbolometer imaging effectiveness is compromised by the 
relationship between PEL pulse width (tpuise= 0.02 ps), pulse period (fpRF’* = 133 ms), and 
the time constant of the pixel absorptive membrane (x = 10 ms). This time constant 
determines the rise time with which the focal plane array responds to incident photons. 
Within a given laser pulse, then, the incident power that is absorbed at any moment by 
the FPA is given by the expression: 

P.,.=P„A-e-"') ( 24 ) 

where t is measured from the start of an individual pulse packet and achieves its 
maximum at tpuise, the total (temporal) pulse width. Because the pulse period is a full 
order of magnitude longer than the thermal time constant, the effects of each pulse are 
essentially identical and independent—that is, sufficient time elapses between pulses that 
total energy absorbed by the FPA pixels decays to negligible levels before the beginning 
of each successive pulse. Maximum laser power absorption occurs at the end of the pulse 
(at t=20 ps); for the laser and microbolometer specifications that have been given above, 
this yields a maximum absorbed power of 12.0 W. However, this is only 0.2% of the 
laser’s peak output power and occurs for an exceptionally brief period of time—for a 
time period that is significantly shorter than the time constant of the detection system. 
The net effect of this short pulse width and slow PRF are evident when one calculates the 
average power which is actually absorbed by the microbolometer over extended periods 
of time (over one or more multiples of one pulse repetition interval): 

= f,,,['’'"P„A^-e-'’)dt = f,„i0.\2mJ) = 0.9mW (25) 
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It is readily seen that—despite its impressively high peak output power—the pulse 
eharaeteristies of this partieular PEL provides an average effeetive power that is no better 
than the QCLs under investigation. Further, as diseussed above, this laser souree also 
earries with it a slow PRF that limits the effeetive video frame eapture rate to only 
7.5 Hz. 

The partieular ease eited here is not a unique one; most other U.S. FEE systems 
have beam eharaeteristies whieh are similarly ineompatible with mierobolometers for 
terahertz imaging applieations. Again, this ineompatibility argues strongly for the need 
for relatively small FEE deviees whieh ean be operated anywhere within the 0.1-10 THz 
range (higher frequeneies, espeeially those eorresponding to atmospherie absorption 
“windows,” being preferable), at high average powers (on the order of a hundred 
milliwatts or more), with relatively high pulse repetition frequeneies (tens or hundreds of 
kilohertz). Sueh a laser souree would allow full exploitation of mierobolometer deteetors 
as an imaging eamera; the two systems together eould likely image metals eoneealed by 
thieker layers of materials already investigated (plastie, rubber, eloth, paper) and possibly 
allow imaging through “new” materials whieh have heretofore been opaque (i.e., 
eeramies, eardboard, ete). 

C. COMPACTNESS 

Mueh of the novelty and allure of this researeh derives from the faet that the 
terahertz deteetor is remarkably portable, is entirely free of eooling requirements, and is 
not dependent upon elaborate optieal elements. Unfortunately, the same ean 
emphatieally not be elaimed for the QCE souree employed in these experiments. 
Although the laser itself is extraordinarily small, the logistieal infrastrueture required to 
support its operation is, for all praetieal purposes, prohibitively eumbersome. Not only 
does the laser require mounting in a stable laboratory environment, its eryogenie 
demands (eooling to nearly liquid helium temperatures) are likewise stringent (adding to 
the system’s already-large spatial footprint). Finally, beeause the QCE output beam 
exhibits sueh a broad divergenee angle, multiple foeusing elements are required to 
eonverge the radiation and produee a eollimated beam for the eamera. (The eryostat 
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window, the presence of which is an unavoidable consequence of inflexible cooling 
demands, adds an additional optical element to the system.) Not until this multi-level 
deficiency in compactness is resolved will a portable system be viable. Fortunately, the 
scientific and engineering communities are actively addressing this shortfall. Indeed, one 
such investigation is well underway at the Naval Postgraduate School. In a novel 
approach to the problem of terahertz source compactness, NFS’s Prof. Andres Larraza is 
currently designing and fabricating MEMS-based, terahertz-tuned micro-magnetrons 
[60]-[61]. Conventional magnetrons (such as those used to generate microwaves) have 
physical dimensions which are excessively large for producing terahertz radiation; typical 
resonant frequencies for magnetrons are roughly 200 MHz. Scaling issues can be 
eliminated by adopting a smaller magnetron structure through MEMS fabrication 
techniques (thus securing a much smaller radius of curvature for electron acceleration); 
however, other technical issues remain. Although maximum magnetron frequency scales 
linearly with applied magnetic field according to the Earmour formula. 


/.= 
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there are practical limitations to the magnitude of magnetic field that can be sustained 
with such a device. In the case of conventional magnetrons, an extraordinarily powerful 
field (approximately 25 T) would be required to generate radiation at 1 THz; application 
of such a field would almost certainly result in an immediate structural failing of the 
magnetron. However, by employing a so-called “reverse magnetron” design (in which 
the polarity of the magnetron elements are reversed—such that the cathode constitutes the 
outer ring; the anode, the central cylinder—a terahertz-resonant cavity can indeed be 
technically feasible [60]. Successful realization of this source technology, if combined 
with the microbolometer detector reported herein, could result in a terahertz imaging 
system that is significantly more compact than any of the real-time systems currently in 
use. 
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D. TERAHERTZ-TUNED MICROBOLOMETER DESIGN 


Having stipulated these eaveats of this system’s technieal limitations, it is 
likewise important to reeognize that nearly all of the teehnieal weaknesses/issues 
assoeiated with the imaging “system” (i.e., low power, fringing effeets arising from 
optieal diffraction, cooling requirements, bulkiness) are liabilities tied to the QCL 
source—not the sensor. To the extent that this research has been primarily concerned 
with the detection mechanism of a potential imaging system, these limitations are of only 
secondary importance to the investigation at hand. Nevertheless, there exist one notable 
(albeit highly work-intensive) opportunity to improve upon microbolometer performance, 
thus optimizing its sensitivity to terahertz radiation: re-engineering of the MEMS design. 

There exist two principal reasons that conventional microbolometer technology is 
not well-suited for imaging at terahertz frequencies. The first is one of physical scale. 
Most IR-tuned microbolometer arrays are comprised of pixels with a 50-pm pitch (that is, 
they have square dimensions of 50x50 pm ). Given that terahertz wavelengths are on the 
order of 100 pm (~ 3 THz), it is plain that pixel pitch should be increased to fully 
accommodate the wavelength disparity between IR and terahertz radiation. This size 
scales linearly with the wavelength ratio, such that a k=100 pm (“terahertz-tuned”) pixel 
should have pitch distance that roughly 5-10 times that of a k = 8-12 pm (“IR-tuned”) 
pixel. This increase would necessitate increasing the pixel length/width dimensions from 
50 pm to 250-500 pm. Furthermore, care must be taken in re-engineering the depth of 
the array for such wavelengths. This is because most microbolometer pixel designs 
tightly engineer the separation distance (typically to one-half the design wavelength) 
between absorbing and substrate layers of the pixel so as to induce a “cavity” effect at IR 
wavelengths (thus optimizing pixel absorption of the incident radiation). However, in the 
case of a terahertz-tuned pixel, it is not feasible to engineer a cavity gap on the order of 
40 pm; therefore, one must use other criteria for determining the optimal depth of the 
sensor channel. 

The second—and more important—element of microbolometer design ill- 
matched to the terahertz regime is that of the pixel membrane. Due to the high TCR of 
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vanadium oxide (VOx) and strong absorption of infrared radiation by silieon nitride 
(Si 3 N 4 ), a eomposite film of these two materials is used as the default membrane material 
used for IR-tuned mierobolometers [55]. And, as this researeh has eertainly borne out, 
the membrane remains at least partially absorptive in the terahertz regime. However, 
despite the faet that Si 3 N 4 absorption has not been well-eharaeterized at terahertz 
frequeneies, eomputational modeling performed by Hu et al indieates that emissivity of 
the pixel membrane at 4.3 THz may be as low as 4% [36]. In light of the imaging results 
presented within this dissertation, this emissivity value seems suspieiously low (such low 
absorption would almost certainly be incapable of supporting the high-contrast imaging 
reported herein). However, it does support the premise that microbolometer operation 
could be significantly improved by adopting a different material for use in the membrane 
film. Identifying a candidate material that exhibits strong terahertz absorption is, in itself, 
a rather complicated prospect; despite the rapid pace of current research into terahertz- 
based technology, published papers on this area of material characterization are relatively 
few in number. A large proportion of the research that has been conducted into this field 
has been performed by the Submillimeter-Wave Technology Laboratory (STL) at the 
University of Massachusetts-Lowell, which specializes in investigations of 0.1-5.0 THz 
[62]-[65]. The STL has measured high terahertz absorption coefficients for Rexolite™ (a 
polystyrene plastic exhibiting 60% absorption at 6.3 THz [63]), hot-pressed hexagonal 
boron nitride (a polycrystalline material, used frequently as an insulating film for metal- 
insulator-superconductor (MIS) structures, with 90% absorption at 3.5 THz [64]), and 
various dichroic fdters, produced using photolithographic techniques, that are currently 
in use as spectral bandpass filters for the terahertz regime and exhibit nearly 100% 
absorption above 450 GHz [65]. Some insight into the problem of material selection can 
also be gleaned from the imaging data presented in this dissertation. For instance, insofar 
as this research has clearly demonstrated that boPET film (such as that shown in Fig. 28) 
is capable of almost complete attenuation of terahertz radiation, Mylar®-type films could 
potentially serve this pixel function as well. 

Unfortunately, even materials demonstrating near-perfect terahertz absorption are 
not necessarily amenable to microbolometer pixel design. While flexibility in choice of 
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material and thickness for MEMS fabrication is continually improving, the current state 
of MEMS production techniques is such that only a relatively modest range of materials 
are available for fabrication of such devices. This restriction serves to eliminate more 
exotic materials from pixel membrane candidates. Still, the exceptionally strong 
attenuation of terahertz radiation by nearly all metal types suggests that basic metallic 
materials (e.g., gold, which has widespread and well-established use in MEMS devices, 
and can be applied in a variety of thicknesses) may ultimately serve as a viable choice of 
pixel membrane medium. (The only immediately apparent disadvantage of this approach 
is the fact that metals have rather low TCRs (i.e., on the order of 0.002 /K, versus 
-0.023 /K for VOx) 

E. RECENT IMPROVEMENTS IN MICROBOLOMETER-BASED 

DETECTOR TECHNOLOGY 

Very recently. Professor Karunasiri’s Sensors Research Group has purchased a 
Thermoteknix Miricle 110 KS infrared (microbolometer) camera for further improvement 
in terahertz detection capability. Notable features of this system (as compared to the 
Infrared Solutions IR-160 camera) include nearly a six-fold improvement in total pixel 
count (384x288 array), a smaller pixel pitch (35x35 pm versus 50x50 pm ), and a 100% 
improvement in NETD for passive, room-temperature imaging in the infrared regime 
(50 mK vs 100 mK). Additionally, the camera is capable of faster data collection (7 ms 
response time and 60-Hz frame rate, in contrast to 12 ms and 30-Hz for the IR-160, 
respectively). The pixel membrane material used in the 100 KS is not publicly known (it 
is considered company-sensitive information); however, it has a design absorptive range 
(7-12 pm) that is slightly more shorter than the 8-14 pm wavelength range that 
characterizes the IR-160. It is also capable of automatically outputting images in the 
form of a colorized intensity map, which may provide better contrast between (and, thus, 
differentiation of) multiple objects that are grouped, overlapped, or concealed within each 
other. 

By all indications, this new camera system should support a measurable 
improvement in the SRG’s terahertz imaging capability. Indeed, preliminary experiments 
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confirm that the focal plane array is markedly more sensitive to the QCL’s 3.6 THz 
illumination than its predecessor. Direct terahertz illumination of the focal plane array, in 
the absence of the Tsurupica focusing lens, produces a noticeably brighter signal pattern 
than was found for the previous detector (Fig. 33). At the time of this writing, the most 
challenging engineering obstacle remains one of optics—successfully integrating the 
Tsurupica focusing substrate into the camera to produce a seamless detection system. For 
optimal performance, the lens must be integrated into the camera system with a lens 
assembly that (1) mates well with the external architecture of the camera device, and (2) 
allows sufficient travel of the Tsurupica element for manual focusing purposes. This 
integration has been accomplished with success in the IR-160 camera, but is still 
underway in the Thermoteknix system. 



Fig. 33. FPA output from Thermoteknix Miricle 1 lOK thermal imager, under illumination by 
3.6 TFlz QCL output beam. The frame shown here makes use of the camera’s contour plotting 
feature, in which varying levels of signal intensity are represented by a range of colors from blue 
(high intensity) to red (low intensity). 


The SRG is currently working on parallel initiatives to (1) ineorporate the 
aforementioned 20-mm Tsurupica lens into the eamera and (2) commissioning production 
of a second, custom-designed Tsurupiea lens (with longer focal length) from Microteeh 
Instruments, Inc. Once eompleted and fully investigated using experimental methods, 
these improvement efforts are expeeted to yield a noticeable enhancement in image 
fidelity. 
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VIII. CONCLUSIONS 


The research presented in this dissertation represents an early but critical step 
toward the realization of a highly robust terahertz imaging system. With the important 
exception of parallel, simultaneous, and completely independent investigations by Prof 
Hu’s MIT group, it also represents the only research into the use of microbolometers as a 
terahertz detection mechanism. Two features of this imaging approach—real-time signal 
detection and uncooled operation—constitute major advantages over the 
currently-adopted techniques that have achieved relatively broad use throughout 
academia and industry. As such, microbolometer-based terahertz imaging is a research 
niche certainly worthy of further examination. In the author’s opinion, further work in 
this area should be pursued aggressively—as should opportunities for academic 
collaboration with other researchers. 

Imaging data presented here confirm the expectation that the strong attenuation of 
terahertz radiation by metallic objects and weak attenuation by nonmetallic objects can be 
exploited to address many important homeland security needs vis-a-vis imaging 
technology. Moreover, by examining a large range of nonmetallic materials, this 
research has more definitively identified the relative ease with which those materials can 
be imaged. To wit, it has been shown that plastics, rubbers, thin paper, and some fabric 
materials exhibit limited terahertz attenuation and are therefore ideally suited for 
terahertz imaging. In contrast, ceramics, foams, thick paper and heavier fabrics prove 
difficult media for terahertz beam transmission at low output powers. It is further 
noteworthy that surprising analytical results have been obtained while imaging domestic 
and foreign currencies (in particular, the heretofore unreported cross-stitch features 
embedded within GBP security threads). These suggest that, in addition to medical and 
homeland security applications, terahertz radiation may hold promise as an anti¬ 
counterfeiting technology. 

Ultimately, a more versatile and portable laser source will be a necessary element 
of any successfully-integrated terahertz imaging system. Given the current state of 
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research into this area, however, more cumbersome sources (such as the QCL reported 
here) will have to suffice for some time in support of basic terahertz imaging research. In 
light of this, the above-detailed evaluation and refinement of techniques for optimizing 
both output power and quality of the NFS QCL beam have considerable value. It is the 
author’s hope that they will serve as a useful operator’s guide for optimizing QCL output 
in future imaging experiments. 

Perhaps one of the most important results from the above research is the clear 
revelation that relatively low laser power is more than sufficient for basic terahertz 
imaging applications. The results obtained using milliwatt-scale laser power (which is, 
incidentally, the lowest average power ever reported in microbolometer-based detection 
experiments; roughly one-fiftieth of the power used in comparable studies by Hu et al 
[36]), holds important implications for similar studies in the future. It confirms that (I) 
high terahertz power is not a prerequisite for imaging success, and, by logical extension, 
(2) given the (eventual) luxury of higher laser powers, much better imaging resolution 
and fidelity should be achievable than was previously expected. 

It would appear, in light of this and many other ongoing experimental 
investigations of terahertz imaging and spectroscopy, that the whimsically-titled 
“terahertz gap” may soon be populated with an array of realized applications. It has been 
soundly demonstrated that many of these applications can be effectively met through the 
use of microbolometer technology. From straightforward imaging of concealed objects, 
to the novel discovery that terahertz radiation can be used to detect heretofore “hidden” 
currency features, to the revelation that low-power laser sources are sufficient for 
realizing each of these important functions, the research presented in the above pages 
confirms that microbolometer-based detection of terahertz radiation holds great scientific 
promise. And, further, that it can provide a viable pathway towards further strengthening 
our national and personal security infrastructure. 
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